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Abstract. We study a spin-i-particle moving on a one dimensional lattice subject to disorder 
induced by a random, space-dependent quantum coin. The discrete time evolution is given by 
a family of random unitary quantum walk operators, where the shift operation is assumed to be 
deterministic. Each coin is an independent identically distributed random variable with values in 
the group of two dimensional unitary matrices. We derive sufficient conditions on the probability 
distribution of the coins such that the system exhibits dynamical localization. Put differently, the 
tunneling probability between two lattice sites decays rapidly for almost all choices of random coins 
and after arbitrary many time steps with increasing distance. Our findings imply that this effect 
takes place if the coin is chosen at random from the Haar measure, or some measure continuous 
with respect to it, but also for a class of discrete probability measures which support consists of two 
coins, one of them being the Hadamard coin. 
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1. Introduction 

Classical random walks are of importance for the field of randomized algorithms |MR95| . e.g. for 
search algorithms, connectivity and satisfiability problems. The generalization of a classical random 
walk to the quantum world, called a quantum walk, is in its simplest form given by a spin- ^-particle 
(qubit) moving on a line, and space as well as time are discrete parameters, see |Kem03| lABN^Ol] 
for reviews. 

Hence the Hilbert space of the particle is given hy n = £2(2) ® C^, where ^2(2) denotes the 
Hilbert space of square summable sequences over Z. Each step in the time evolution is described 
by the same unitary operator, called the walk operator. It is defined as the product of two unitary 
operators, the first one being a shift operation which acts on the position degree of freedom of the 
particle. The shift operator S moves the particle one position to the left or to the right, depending 
on its internal degree of freedom, 

S{5x e±) = 5x±i e± 

where e+ and e_ label an orthonormal basis of C^, e.g. spin up and down states, and Sx denotes the 
state of a particle which is localized at position x on the lattice. Mathematically speaking, Sx is just 
the element (. . . , 0, 0, 1, 0, 0, . . . ) of the canonical basis of ^2(^) which has only one non-zero element 
corresponding to the point rr G Z. The second operation is a unitary transformation acting locally 
on the internal degree of freedom, defined by the direct sum of elements of the unitary group in two 
dimensions, 

U = ^Ux, Uxe U{2) . 

x^TL 

We speak of the \Jx as quantum coins or just coins, and use the term coin operator for U. It is 
block-diagonal with respect to the identification of ^2(^) ® with ©ie^C^. If the coin action does 
not depend on the position of the particle in the lattice, Ux = V for all x € Z, then U has the simple 
form [/ = The total walk operator is now defined to be the product 

W = U-S. 

Of course, interesting effects can only take place if the coin creates a superposition in the basis 
elements e± of C^, so that the particle is shifted both ways. For translation invariant walks Fourier 
methods have been a fruitful method to determine the position distribution for finite times as well 
as asymptotically |AVWW10| . 

It is well known that a quantum walk can exhibit a propagation speed which is quadratically faster 
compared to a classical random walk |Amb03| . Hence, it is natural to ask for their use in the theory of 
computational algorithms. And indeed, there are several proposals for quantum algorithms making 
use of the concept of quantum walks |Kem051 IAmb03| . which improve on their classical counterparts. 
Another feature of quantum walks is that experimental realizations are challenging but feasible with 
current technology |KFC+n9[ [SMS+n9| . 

Such implementations of quantum walks, however, involve the control of many experimental pa- 
rameters, e.g. laser beams or microwaves. Since this can only be done with finite accuracy, noise will 
be introduced to the system. There are several extremal cases of noise, for instance we could assume 
that the control parameters vary homogeneously in space and on small time scales compared to the 
execution time of the quantum walk. This means, that the time evolution changes over time, i.e. 
the coin operation becomes time dependent. This kind of noise has been studied in |AVWWTO] and 
it was shown that this leads to diffusive behavior of the quantum walk, i.e. the walk behaves like a 
classical random walk. 

Here, we consider the scenario where the control parameters vary on large timescales compared 
to the execution time of the experiment, but spatial disorder breaks translation invariance. This 
situation is modeled by a random, space dependent coin. These fluctuations resemble a random 
potential, but are still described by a unitary operator. The question we address is whether the 
performance of quantum walks is seriously affected by such kinds of disorder. From the theory of 
disordered crystals we would expect to observe localization phenomena similar to those described 
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by Anderson in |And58| . There, imperfect crystals are also modeled by a random, space dependent 
potential. The strength of the varying potential is specified at each point by some probability 
distribution, which is usually assumed to be independently and identically distributed, i.e. it is 
the same at each point in space. The Hamiltonian describing such a system is then a sum of the 
lattice Laplacian and a random, multiplicative potential. Anderson argued that due to interference 
effects, electrons in disordered crystals can exhibit a strong tendency to localize in finite regions up 
to exponential tails. In other words, the tunneling probability between two lattice sites decays for 
all times with increasing distance faster than any polynomial. This phenomenon is called dynamical 
localization |Kir07| . 

We show that the statement still holds for the case of a spin-^ disordered quantum walk in 
one lattice dimension, where the coin is assumed to vary probabilistically over space, but remains 
fixed for all times, under mild conditions on the coin distribution. The class of measures for which 
these conditions are verified include distributions of coins continuous to the Haar measure, e.g. a 
gaussian distribution with small variance and centered around some coin of interest, but also a class 
of probability measures supported on two unitaries, one of them being the Hadamard matrix. 

Related models have been studied by Joye and Merkli |JM10| and Hamza et. al. |HJS09| . which 
can be seen as a study of disordered quantum walks with identically independent distributed random 
phases, chosen with respect to an absolutely continuous probability measure. Hamza and co-authors 
considered diagonal disorder, and find that localization effects occur only if the shift is assumed to be 
imperfect. Joye and Merkli generalized these results to the case where the coin is given by a product 
of a fixed unitary with a diagonal one with independent random phases and observe localization 
effects also for the deterministic shift. Linden and coworkers |LS09| considered a case where the 
perturbation of the coin operator is periodic in space, and Konno and coworkers |Kon09al IKon09b| 
as well as Shikano and Katsura |SK10| derived limit theorems for a special kind of disorder. Compared 
to previous work, we consider all possible coin distributions and derive sufficient conditions such that 
they lead to dynamical localization. 

Our proof incorporates ideas from the study of localization effects in random unitaries |HJS09| 
IBHJ03I Joy04| , and uses some general results about Borel measures on the unit circle (see the Book 



by Cima, Matheson and Ross |CMR06] for a really nice treatment of this theory) . We rely heavily 
on the theory of products of real- valued random matrices developed mainly by Fiirstenberg |Fiir63| . 
which we use as a starting point to derive limit theorems of products of complex matrices having 
determinant one. As a note for specialists concerning localization questions, let us mention that the 
proof is based on a multiscale-analysis-scheme, and we first derive the necessary statements, i.e. a 
Thouless-like formula and a Wegner estimate. 

The paper is organized as follows. We first define in a mathematically rigorous manner the notion 
of a disordered quantum walk and state our main theorem. The proof is split into different lemma 
and is given in section HI Some technicalities as well as the verifications of our conditions for our 
examples of coin distributions are shifted into the appendix. 

2. General Setting 

A disordered quantum walk is a quantum walk where the time evolution is defined by an element 
of a family of random walk operators W^j with space dependent random coins. That is, the shift 
part S of the walk operator is assumed to be undisturbed, but the action of the coin operator at 
each lattice site x E Z is now given by a random variable taking values in the group U{2) of two 
dimensional unitary matrices. 

We will call the distribution of this random variable the single site distribution and will denote it by 
jjLx- From now on, we call the set 

supp(^,) := {M e GL(C, k) : ^ix{Be{M)) > OVe > 0}, 

the support of fx^- Here, Be{M) denotes the open sphere of radius e around M given by the operator 
norm. We require that the coin operations are independent and identically distributed random 
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variables at each lattice site, that is = 1^- Correspondingly, the joint distribution of coin operations 
for a finite collection of L lattice sites is given by the product of L copies of the single site distribution 
fi. Since this is true for all finite collections of lattice sites, there exists a unique joint distribution 

fioo on the infinite product probability space 

^ = Xxez ■ 

We will use the symbol IE {X) to denote the expectation value of some random variable X, and F [E) 
to denote the probability of some event both taken with respect to /ioo- The coin operator 
a; S f7 on £2(2) ® is a random variable on $7 and defined by the direct sum 

= t^(...,tj^_i,tj^,a;j;+i,...) = ^3 . 

The one-dimensional, spin-^ disordered quantum walk W^j is a random variable on with values in 
the unitary group of ^2(2) ® given by the product 

W^ = U^-S. 

Note that the coins are drawn in advance and kept constant during all time steps. 

The coins of a quantum walk play the role of a potential, i.e. they change the probability that a 
particle is moving to the right or to the left, therefore determining the transmission probability. If all 
coins are the same, i.e. the system is translation invariant, the generalized eigenfunctions are described 
by Bloch waves being infinitely extended over the whole lattice. Thus, they are of course not square 
integrable and hence no elements of the Hilbert space. Because of that, they are called generalized 
eigenfunctions. The associated spectrum of the Walk operator does not consist of separated points, 
but is rather constituted of energy bands. It follows that the walk operator has absolutely continuous 
spectrum and the system exhibits ballistic scaling. See |AVWWTO] for a thorough treatment of this 
case. 

As the coins vary in space, the system is no longer translational invariant, and the particle "sees" 
different "potential values". Due to that, we would expect that the spreading of initially localized 
wave packets will be slowed down. For the case of hamiltonian systems, we even have that the 
tunneling probability between two lattice sites is suppressed in their distance. This phenomenon is 
called dynamical localization. The corresponding definition for quantum walks is as follows. 

Definition 2.1. Let be a one-dimensional, spin-^ disordered quantum walk. W^j is said to exhibit 
dynamical localization, if there exists a function £ : N — t- 1R,+ such that C{n) goes to zero for n — t- 00 
faster than any polynomial in n and we have that 

^(sn^\{5y®(l),W%®il>)\] < C{\x-y\). 

Here, 5x and 5y denote again localized (position) states, whereas (j) and i{j refer to some arbitrary 
internal spin states. The function C is called the localization length. 

Dynamical localization is closely related to the spectral properties of W^j. As a quantum walk is 
described by a unitary operator, its spectrum a(W;^) is a subset of the unit circle T = {z € C, \z\ = 
1}. We will call its elements frequencies or quasi-energies. The disordered quantum walk is said 
to exhibit spectral localization if it has only pure point spectrum, meaning that each of its spectral 
measures is a linear combination of Dirac point measures. Due to a version of the well-known RAGE 
theorem adapted for unitary operators (theorem ID.3P dynamical localization also implies spectral 
localization. It follows that Wi^ has a complete set of square integrable eigenfunctions. In the case 
of disorder, we would expect even a stronger decay behavior of eigenfunctions. 

Definition 2.2. A disordered quantum walk W^j is said to exhibit strong spectral localization, if its 
eigenfunctions decay faster than any polynomial, implying that the spectrum is of pure point nature. 
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3. Results 

To state our results concerning disordered quantum walks, we first introduce the following family 
of mappings from the general linear group in two dimensions GL(C, 2) into itself 

a b \ 1 ^ 



c d I a \ —b z 



c 



where of course a 7^ is assumed. In the special case when the mappings are applied to unitary 
matrices, it can easily be seen that their image consists only of elements of the group SL^ of complex 
matrices with determinant of modulus one, 

SLt = {M G Mat(C,2) : |det(M)| = 1}. 

Here, Mat(C,2) denotes the space of two by two complex matrices. The matrices Tz{U) G SLt, 
where U £ U{2) is in the support of fi, are called transfer matrices. Their importance stems from 
the fact that the spectral and dynamical behavior of disordered quantum walks crucially depends on 
the group of matrices generated by the set of transfer matrices. 

Theorem 3.1. Let he a disordered quantum walk characterized by its single site distribution fi 
with support supp{fi) C U{2). Suppose that for almost all elements of the unit circle G T the group 
(/ie) C SLt generated by its associated set of transfer matrices Tg{U), U £ supp{(j) is 

(1) non- compact, 

(2) contains no reducible subgroup of finite index and 

(3) the expectation IE ^| |t0([/)| is finite for some > 0. 

Then W^j exhibits dynamical as well as strong spectral localization. Here, the terminus "almost all" 
refers to the Lebesgue measure on the unit circle. 

Remark 3.2. // the group {^q) C SLt happens to be compact for some element 9 of the unit circle, 
then 9 would be an element of the absolutely continuous spectrum ofW^^. This is exactly what happens 
in the case of translation invariant quantum walks. 

This result implies dynamical localization for many physical situations, i.e. for a setup where the 
application of a fixed coin is desired but it cannot be circumvented that with some small probability 
quantum coins close to the target one are also applied. 

Corollary 3.3. Let be a disordered quantum walk such that the single site distribution possesses 
a positive density with respect to the Haar measure onU{2). Then W^j exhibits dynamical as well as 
strong spectral localization. 

In fact, it is enough if the single site distribution consists partly of some absolutely continuous 
measure. The verification of the assumptions of theorem 13. II for this case is carried out in Appendix 
[B1 Although continuous distributions can be argued to model many experimental imperfections, 
the situation where the support of the single site distribution includes only a finite number of uni- 
tary matrices is of independent interest. Indeed, we may ask what happens if we disturb the usual 
Hadamard walk, as realized by |KFC+n9[ [SMS+n9| . by just one single coin chosen with some arbi- 



trary, but non-zero probability. Then, we have the following result, as shown in Appendix [Q 

Corollary 3.4. Let be a disordered quantum walk such that the support of its single site distri- 
bution equals the set 

' J_ / 1 -1 \ (a h 
V2 U I ) '\-h a 

with complex numbers a and h fulfilling a 7^ 0, |ap + \h\^ = 1, and \a\ < \'3h\, where 3h denotes the 
imaginary part of h. Then W^^ exhibits dynamical as well as strong spectral localization. 

In the case that a = 0, there is a non-zero probability that "flips" occur. These are coins which 
cause reflections of the walking particle at the corresponding sites, see section 14.41 Dynamical 
localization follows easily in this case from lemma 14.81 We continue with the proof of theorem 13.11 
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4. Proof 

4.1. Outline of the proof. At first, we reduce the problem to the case of finite lattices and show 
how the expectation value of interest can be bounded by an expression involving the resolvent of the 
quantum walk operator, see equations ([I]), (l6|) in section 1^2] The next sections. 14.31 and 14.41 study 
the decomposition of the resolvent into expressions mainly depending on transfer matrices, see in 
particular lemma 14.91 In section 14.41 the properties of transfer matrices are examined, and we use 
the theory of Fiirstenberg, (appendix |A]), about products of random matrices to get a first estimate 
of the decay properties of the tunneling probability for a fixed length scale (proposition I4.l4|l . We 
then proceed by examining the properties of the density of states, both for the finite and infinite 
lattice case (section 14. 5p . In particular, we prove an analog of the Thouless formula for quantum 
walks (see equation pS]) ). The Thouless formula is then used to prove the Holder continuity of the 
integrated density of states (see proposition I4.2ip based again on general properties of products of 
random matrices (appendix |A]). This in turn allows us to prove an upper bound on the probability 
that the tunneling probability for some fixed quasi-energy is high for two independent regions at the 
same time, i.e. we prove a Wegner-type bound (proposition 14. 23) and equation ([28]) ). The initial scale 
estimate and the Wegner bound are then combined using the multiscale-analysis technique, similar 
to the ones in |GKOH IKle08| , in order to get an upper bound on the tunneling probability valid for a 
sequence of increasing distances, see lemma 14.241 The final section, 14.91 combines these results and 
finally provides an upper bound on the decay properties of the tunneling probability over arbitrary 
distances, as well as bounds on the decay of eigenf unctions. 



4.2. Restrictions to finite volumes and Cauchy transform. In this section we define a finite 
unitary restriction Wuj{N) of the Walk operator on the lattice sites —N to N. The idea is to 
impose refiective boundary conditions on the lattice sites —{N + 1) and A'^ + 1 which is equivalent 
to changing the coins Ui at these two lattice sites to the flip operation, e.g. the Pauli X matrix 



X 



1 

1 



multiplied by a phase factor e*'' ' . Then we restrict the operator to the lattice sites — to A^, but 
in order to obtain an unitary operator, we, roughly speaking, have to include half of the operator on 
the neighboring sites — (A^ + 1) and A^ + 1 where we changed the coin to be the flip operation, see 
also (d]). 
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To be more precise we define the restriction W^{N) : C^(^+i) ^ (C^iN+'^) of exphcitly via its 
matrix elements 

(2) {WUN))i,j = {W^h,j, -2N-l<i,j<2N + 2, 

where is the walk operator with the coins f7±(Ar+i) changed to the Pauli X matrix multiplied by 

a phase factor e^v^'^ . Note that this choice already ensures unitarity of Wuj{N) since all rows and 
columns are normalized and mutually orthogonal. 

Since the only transport between neighboring lattice sites is caused by the shift S the walk operator 
Wui as well as for finite t G N are by construction band matrices of finite width (of course 
depending on t). This implies a finite velocity for any initially localized particle, meaning that we 
can provide a t-dependent upper bound on the region where the particle could be detected after t 
time steps with non-zero probability. Hence, we can express the probability to detect an initially 
at position x localized particle at position y after applying a fixed number t of time steps by finite 
restrictions of quantum walk operators, 

(3) \{6y(S>cP,W^6x(S>^p)\ = lim \ {6y (g) ,WUNYSx (g> ^P)\ . 

Here (f) and ip denote again some arbitrary internal spin states. But clearly, it is enough to consider 
the case (j) = Cj, ip = ej, i,j E {1,2} where ei, 62 is the standard basis set of C^. We make another 
simplification and consider the time behavior of Wuj{N) first just for open subsets of quasi-energies, 
i.e. open arcs of T. That is, we prove the decay of quantities of the form 

I {6y <S) Ci ,W^{Nyx{Is)6x (S) ej)\ 

where denotes an open arc of the unit circle centered around some 9 G T with elements whose 
phases differ from the phase of 6 only up to (5 > 0, and xi^s) is the associated (eigen-)projector. Once 
we proved the decay for such an open arc, the result for the whole unit circle follows by compactness. 

Let us denote the spectral measure of W^(A^) associated to the vectors 6x Cj and 5y Cj by 
Puj^N- omit here the dependence on ei,ej, since all our arguments are independent from them. 
Note that since W^(iV) is a unitary operator acting on a finite dimensional Hilbert space, its spectral 
measures are singular measures, i.e. only supported on points. By the spectral theorem and since 
Wuj{N) is unitary, p^'^j^f is a signed measure on the unit circle T and 

sup II {6y (g) e^ , W^{NYx{Is)6x «> ej )\\ = sup 
ieiN tev 

where we denoted in a slight abuse of notation also the characteristic function of the set Is by xi^s)- 
Since \6^\ = 1 for 9 gT we can bound the expression on the left hand site using the triangle inequality 
by the measure of the open arc Is, 

sup|(5,®e,,VF,(iV)*x(W.®e,-)| < / \p:,'l{d9)\x{Is) = Ip^'^m ■ 

In order to estimate the value associated by p^'jy to the open arc I5 we introduce the Cauchy 
transform Kp^'^j^, defined as a function on the open unit disk B = {z G C : |z| < 1} by the formula 

kp:%{z) = [ p:%id9)-^. 

Jt I - 9z 

Note that Kpl]\{z) is an analytic function on the open unit disk. It may become infinite for some 
values of z as z approaches the unit circle, and these exceptional values are precisely the points where 
the spectral measure is supported and hence are equal to the eigenvalues of the walk operatoiQ. A 
theorem of Smirnov assures that except from these points (which have Lebesgue measure zero), the 
limits \\m^^i- K p^'^j^{r9), 9 GT exist. See |CMR06| for a formal statement of this theorem as well 
as for a proof. The reference is also well suited for a general introduction into the topic of Cauchy 



/ p:^'%{d9)xii5)0' 



Note that we are dealing here with unitary operators on finite dimensional Hilbert spaces and hence the absolutely 
continuous part of every spectral measure is identically zero. 



8 



ANDRE AHLBRECHT, VOLKHER B. SCHOLZ, ALBERT H. WERNER 



integral transforms and their properties. Shortening the notation, we denote the corresponding, 
ahnost everywhere weh defined function on the unit circle again by Kp^'^j^. 

A theorem by Poltoratski |Pol96| . which can be thought of as a refinement of Boole's well known 
theorem, states that we can express the action of a singular measure on the unit circle on some 
continuous function / by studying its Cauchy transform, 

(4) I \p:Ud9)m = Jim .K.|^dex{,,.,^,^,y,)|>.}W/W 

where dO is the integration with respect to the normalized Lebesgue measure on the unit circle and 
XA again denotes the indicator function of the set A. Put differently, the measures 

TTKdex{ e&J:\KpZ%{e)\>K} 

converge to the absolute value of the spectral measure p^\{dO) in the weak-*-topology as k goes to 
infinity. Remarkably, this theorem also holds if we restrict ourselves onto some open arc of the unit 
circle, which gives us 

(5) \p::\m = Jim vrK.A|{0G/5 : \Kp:^y^{e)\ > 

Here, \\A\ = dO denotes the Lebesgue measure of the set A. Hence we have to control the growth 
of the Cauchy transform of the spectral measure p^'^j^- 

In the special case when the measure on the unit circle is induced by a unitary operator, the 
Cauchy transform is directly linked to the corresponding matrix elements of the resolvent, 

x,y / \ / x,y 



= {dy (g) ei , WUN){WUN) - z)-^6^ ® ej ) 
= {dy(S)ei,l + z- {W^{N) - ® ej ) 

= z- {5y®ei, {W^{N) - ®ej) , 

for X ^ y. Hence, we have that for x ^ y 

(6) \pl;%\{h) = yim^^^-\\{e^h ■.\{5y®ei, {W^{N) - Oy'S, 0e,)\ > k }| . 

We continue by first deriving an explicit decomposition of the resolvent in terms of the already 
mentioned transfer matrices, see the next sections. In section we study their properties and we 
finish by proving the decay for some open arc Ig in section 14.91 



4.3. Resolvent formula. In the following we derive an explicit expression for the resolvent in terms 
of transfer matrices. Since the results are independent of the size N of the quantum walk W^j^N), 
and for shortening the notation, we mostly omit the dependence on N. We can identify the Hilbert 
space £2(2) ® with ^2(2) via the isomorphism 

4 (g) e_ f2x , 5^(S>e+ ^ f2x+i , G Z , 

where the e± form an orthonormal basis of and the fj denote the standard basis of ^2(2)- We 
want to study the matrix elements of the resolvent with respect to the standard basis vectors 5x0 e±. 
To shorten the notation let us abbreviate the resolvent by 

G. = {W^ - z)-^ , z i a{W^) . 

We denote the Tth column vector of the resolvent by 7' G ^2(^)- In the following we determine 7' 
locally in the sense that knowledge of two adjacent entries 723, and J2X-1 allows us to construct 7' 
iteratively. Assume now we know the entries 723, and ^2x-i ^'^'^ ^ 2x,2x + 1. The shift and coin 
operators acting on ^2(2) are given by 



Sf2x+i 



f2(x-l) ,1 = 

f2{x+l)+l 5^ = 1 
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2x+i 





axf2x + Cxf 2x+l 
bxf2x + dxf2x+l 



, x = y, i 
, x = y, i 



Hence, we can determine the operatoi0 



—z 





a-x-i 








> 







—z 


Cx-l 














bx 


—z 





ax 








dx 





—z 


C-x 














bx+i 


—z 














dx+i 





—z 



V 

and it follows from the equation 

{W^ -z)-^' = fi 
that for / 7^ 2x,2x the column vector 7' satisfies the equations 

l2x-ibx - z-iL + l2x+2ax = 

l2x-ldx - + l2x+2Cx = . 



C-x 

z 



l2x-l 
7L 



It is easy to see that this implies the following relation between the entries T2x-i^ i2x i2x+i^ 
7L+2 

I / det(j7^) 
^2x4-2/ «^ V -bx 

which motivates the following definition. 

Definition 4.1. Let W^, = U^j ■ S be a disordered quantum walk. The transfer matrices Tx{z) ofW^, 
for all sites x with a^; 7^ are defined by 

I / det(U^) 
-bx 



Tx(z) :-- 



C-x 
Z 



x € Z. 



Hence, with the definition 



r 

J- r. 



7L-I 
7L 



and under the condition I 7^ 2x, 2x + 1 we have the relation 

r'L+i = Tx{z)T^^ . 

Remark 4.2. The condition Ox ^ is crucial for the transfer matrices to be well defined. In fact, 
the case Ux = corresponds to a reflection of the walking particle at site x. We analyze this case in 
the next section and show that we can assume without loss of generality Ox ^ for all sites x. 

If the requirement I 7^ 2x, 2x + 1 is not satisfied, the application of Tx{z) to the vector results 
in a vector different from ^x+i- '^^^ following lemma is concerned with an analysis of the sequence 



iLx 



defined as 



(7) 



d,x 
4>2y 



Ty{z)-K..Tx-i{z)-^T^, ,iiy<x 

, if y = 2; 

Ty^i{z)...Tx{z)T[. ,iiy>x 



The columns resp. rows of the matrix are indexed by integers increasing from left to right resp. top to bottom. 
The block in the middle corresponds to basis vectors f2x and f2x+i. 
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Lemma 4.3. Let y,x,l € Z and (plf be as defined in The resolvent and the sequence (j)^ 

coincide, i.e. 

if2x — l<l and 2y — 1 < I respectively 2x > I and 2y > I. 

Proof. The lemma follows directly from definition 14.11 □ 

Remark 4.4. In the case of finite restrictions Wuj{N) the sequences (plf have to fulfill certain 
boundary conditions. In the limit N ^ oo it is clear that the resolvent is a bounded operator and 
therefore it follows from lemma\4^ that for 2x — 1 < I the sequence (j)^ is left square summable, i.e. 



and for 2x > I it is right square summable, i.e. 

n>r 

In the following we call sequences (p± satisfying these requirements left respectively right compatible 
with Gz. 

In fact, since the series ([7]) is a solution of the equation 
(8) {W^ -z)^ = 

for z ^ (t(W;^) it can only be left or right square summable but not both. Now, let <j)-{n) and (j)+(n) 
be left (respectively, right) compatible solutions of ([8]), for instance the ones given by lemma If 
we choose exactly the solutions of lemma 14.31 it is clear that the matrix element (n,m) is given by 
(p-(n) respectively <j)+{n). For general left respectively right compatible solutions (j)± we write for the 
resolvent at m = 2y,2y + 1 



Gz (n,m) 



(3m(l)+{n) , ifn>2y. 



When we consider the cases m = 2y or 2y+l the constants am and have to satisfy the equations 

am{by4>-{2y~l) - Z(j).{2y)) + f3may(l)+{2y+2) = 5m,2y 
amdy<i)-(2y~l) + Pm{Cy(p+{2y+2) - Z(i)+{2y+l)) = I - 5m,1y ■ 

Since (/>+{■«) and (/)^{n) are solutions of ([8]) we can rewrite these equations in the following form 

-ay(j)^{2y+2) ay(j)+(2y+2) \ f Oim \ ^ f ^m,2y 
dy(\)^{2y-\) -dy(t)+{2y-l) ) \ Pm J \ - 5m.,2y 

Using Cramer's rule and the definition 



4>~(2y-l) (t)+(2y-l) 



y ' (t)~{2y+2) (t)+{2y+2) 

we see that the coefficients and /3m are given by 

_ <p+{2y~l) a _ 0-(2y-l) 

— aydet{Ay) ' P'^V ~ aydet{Ay) 
_ <P+{2y+2) o _ 0-(2y+2) • 

"2J/+1 — dydc\.{Ay) ' P2y+l — dydet{Ay) 

We would like to replace Ay by a matrix By where the components of (j)± appearing in By correspond 
to adjacent sites. That would allow us to calculate By from via transfer matrices. Arbitrary 

solutions (j) oi ^ satisfy the relation 

0{2y-l) \ M ^ \( '/'(2?/-i) 
(l>{-iy) ) V T ^ / V '^(2y+2) 
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hence, we have for 

(9) By := ( ^f'-''' 'f'f'-'^ 

the following correspondence 

det(Sy) = ^detiAy) . 
Moreover, we have By^i = Ty(z)By and since det(rj^(2:)) = dya"^ the matrices By satisfy 

|det(53,)| = |det(5y)|, Vy,y'GZ. 
Finally, we get the following expressions for the coefficients am and 



zdct{By) ' ^2j/ — zdct{By) 



_ ay<p+{2y+2) o _ ay<P^{2y+2) ' 

"22;+l — zdydct{By) ' P'iy+^ — zdydct{By) 

The following lemma summarizes the results of this section. 

Lemma 4.5. Let W^j be a disordered quantum walk and (j)-{n) and (j)+{n) left respectively right 
compatible solutions of Then the modulus of the entries of the resolvent (n,m) for m E 

{2y,2y + 1} is given by the formula 

nn^ IP . J _ 1 / \4>~M4>+i^)\ , n<2y 

^ ^ ' ^ " \zdet{B,)\ ' 1 |<A+(")</'-(-)| , n>2y 

with Bk defined by ^ and arbitrary /c G Z. The argument m is given by m = m — 1 for even m and 
rh = m + 1 for odd m. 

4.4. Transfer matrices and properties of the resolvent. In this section we analyze the resolvent 
formula (|10p in further detail and prove some useful facts originating from the structure of the transfer 
matrices T{z). Those matrices induce a family of maps labeled by z € C\{0} from a subset of the 
unitary group U{2) to the group of invertible matrices GL(C, 2) via 

(11) : ( ^ ^) e^(2),a/0 i ( ^ ^ ) G GL(C, 2) . 

In order to make the map well defined, we define the set of two by two unitary matrices with non 
vanishing diagonal entries 

(12) Und := {U G U{2) : f/n / / U22} , 



the complement of Und in U{2) is denoted by Und = U{2)\Und- For each U G Und we have that 
det(rz(C/)) = U22/U11, hence it is clear that Tz(Und) C GL(C,2). The following lemma proves that 
this inclusion is strict and gives an explicit parametrization of Tz{Und)- 

Lemma 4.6. Let Tz and Und be defined as in Ul\) and Then, for z G C\{0} is injective 

on Und OLf^d its image Tz{Und) is given by the set 



(13) [[ ) : rGR+,a,/3,7G[0,2vr) 

Proof. The inverse of is given by 

-I f V w \ ^ f ^ ~^ 



(14) t; , 

' X y J y y w vy — xw 

which only exists for y ^ 0. Since we assumed z ^ this already proves that o Tz{U) = U for all 
U G Und, hence is injective. 

The image of can be inferred from its action on a general unitary 



Clearly, if we define r = |6|/|a|, then r G 1R,4- and \a\ ^ = \/l + r^. The phases a,/3 and 7 are 
arbitrary, because (p and the phases of a and b are arbitrary. □ 
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The parametrization (jl3p immediately gives the following invariance property of the transfer ma- 
trices. 

Corollary 4.7. Denote by P2 the set of vectors x G such that 

X=(^MgP2 <^ |xi| = |x2|, 

then for z G T and T(z) G Tz{Und) we have T2P2 C P2- 
Proof. Without loss of generality we may assume 

X = ^ ^ with |c| = 1 . 

Now, the statement follows from the fact that 



iVl + rV" + re'^c\ = \re'^ + yjl + r2e*(^+^-")c| . 

□ 

Since is injective on IAmd^ a measure /x on ^(2) induces a family of measures /i;^ on Tz{Und) 
via 

If /i is a probability measure on ^(2) the fi^ are probability measures on Tz{Und) if and only if 
fi{Ul\f£)) = 1- In this case, with probability one, Uaj{N) is a direct sum of elements in Und such that 
the transfer matrices r(^) for Wuj{N) are well defined. Consequently, it is valid to assume that the 
resolvent {Wui{N) — z)~^ can be expressed through transfer matrices T{z). 

The case < 1 also leads to dynamical localization of W^(A^). The reason is that coin 

operators Ux, Uy G Und cause reflections of the walking particle at sites x and y, see |LS09| . hence 
a particle starting in between x and y is strictly localized in the flnite region between x and y for all 
times. The following lemma proves dynamical localization for iJ,{Ui^u) < 1. 

Lemma 4.8. If ijl{IAnd) = p < 1 the disordered quantum walk W^j satisfies 

lE{\{6y(E)c|>,WX^^P)\) < {l-p)\-y\+^ 

Proof. We assume without loss of generality x < y. Let x < r £ 1, denote the site such that 
Ur G Und and Ug G Und for all x < s < r. If Ug G Und for all s > x we set r = 00. The scalar 
product satisfies 

\{6y(g)cl),WX^i^)\ = 0, Vy>r,Vt>0 

and the probability that Ur G Und and Us G Und for all x < s < r is exactly (1 — pY~^~^p. Hence, 
the bound of the lemma follows by assuming the worst case, i.e. 

I {5y^(l),W%0^)\ = 1 

for all y with r > y. □ 

The preceding lemma tells us that we may assume that all coins Ux of the walk operator W^^ 
satisfy Ux G Und when proving dynamical localization for W^j. This leads us to a simplification of 
the resolvent formula (|10p for finite walk operators Wi^{N). 

Lemma 4.9. Let he a disordered quantum walk and Wuj{N) a finite restriction of it. Denote by 
the resolvent 0/ W^(A^), then for z G T there exist normalized vectors ^>4- G P2 such that 
the matrix elements {2x-i,2y-j) for i,j G {0, 1} obey 

^ ^ ' ^ ' 2 \ \{^+,Tx-iiz)...Ty{z)<!>^)\ ^ , ifx>y 

The vectors ^± depend in a non-trivial way on N, x, y and z. 
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Proof. For the moment we suppress the dependence on z and denote by cp^ and (j)^ the left respec- 
tively right compatible sequences for . Our choice of boundary conditions implies 

\4>^i-2N-l)\ = \(f)^(-2N)\ and |(/)^(2Af+l)| = |(/)^(2Af+2)| , 

and hence, by Corollary K7\ we get \(j)^(2M-i)\ = \(j)^{2M)\ for aU -N < M < N + 1. Now, in the 
case X < y the statement follows from (fTO|) if we set 



^>_ = —= — and $ 



V2|0^(2x)| V 'P^i^x) J + V2\cl>^i2y)\ \cl)^{2y-l)J ' 

and in the case x > y we have to exchange x and y in this definition. □ 

Remark 4.10. We note in particular that given two normalized vectors (pL, (pn G P2 and a set 
of transfer matrices, then we can always construct some restricted walk operator such that the right 
hand side of (fT5|) constitutes the resolvent of that unitary operator. This follows by choosing the right 
boundary conditions in the construction of the finite restriction outlined in since the mapping of 
unitary coins to transfer matrices can be reversed. 

We also need the following estimate on the difference between products of transfer matrices of 
length n associated to two different elements of the unit circle. 

Lemma 4.11. Let n G N, G T and assume that the transfer matrices Ti, . . . ,T„ satisfy \ai\ > n. 
Then 

F {\\\Tn{e) ■ . . . ■ T,(e)v\\ - \\Tn{e') ■ . . . ■ TAe')v\\\ > r]) < 

ri 

for all normalized vectors v G C^. 

Proof. By Markov's inequality, the inverse triangle inequality and the definition of the operator norm 
it is sufficient to prove 

E (||r„{e) • . . . • T^e) - r„(e') • . . . • T,{e')\\) < C{7i)\e - e'\ . 

Writing the operator difference as a telescope sum and using again the triangle inequality we obtain 
the bound 

WTniff) ■ . . . ■ Tm - Tnie') ■ . . . ■ Tiie')\\ < 

\\Tnie)-Tn{9')\ \ -WTn-m- ...-TmW + 
\\Tni9')\\ ■ \\Tn-m - Tn-li9')\\ • | |T„_2(e) • . . . • Ti (0) 1 1 + 



+ 



\\Tn{e')...T2{e')\ \ ■ \\Ti{e) - Ti{e')\ 
For products of transfer matrices we have the operator norm bound 



2\ A: 



\\Tk{e)...Tm\\<(-] , V^gC, 1^1 = 1 



which follows from the norm inequality ||A • < \\A\\ ■ \ \B\\ and the fact that the singular values 
X± of Ti{e) for 16*1 = 1 satisfy 

A± = ^ < - . 

The assertion of the lemma with C(n) = n-2"~^/K" follows now from the equation ||Tj(6») — Tj{e')|| = 
\9-e'\/\a,\. □ 

We also need the following easy corollary to lemma B. Ill 
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Corollary 4.12. Let n G N, ^ € T and assume that the transfer matrices Ti, . . . , r„ satisfy \ai\ > k. 
Then 

F ( I I {v^,Tn{e) • . . . • T^{e)v^)\ - \ {v, ,Tn{e') T^{e') v^)\\>ri)< -C{n)\9 - e'\ 

V 

for all normalized vectors vi,V2 G C^. 

Proof. After applying Markov's and the inverse triangle inequality we can just use the Cauchy- 
Schwarz inequality for the scalar product. Since the vector vi is normalized we arrive at the same 
situation as in the proof of lemma 14.111 and the whole argument carries through also in this case 
with the same bounds as before. □ 

In theorem 13. II we assumed that the group (/x^) generated by the transfer matrices T{z) G supp(/i) 
is non-compact and possesses no reducible subgroup of finite index, for some z = 9q G T. The reason 
is that we wish to apply results of |Fiir63| to show that with high probability the scalar product 

I {vi,Tn{z) ■ . . . ■ Ti{z)V2)\ 

grows exponentially for all normalized vectors vi and V2- More precisely, we prove the following 
corollary building on |Fiir63| in appendix |Al To simplify notation we set x = 1 and y — 2 = n for 
the remaining part of the section. 

Corollary 4.13. Let fi be a measure on hl{2) and suppose that the group (/i^) generated by the 
induced measure is non-compact and possesses no reducible subgroup of finite index for some 
z G C, then the product 

\\Tn{z) ■ Ti{z)v\\ 

grows exponentially in n almost surely. Furthermore, there exists 7, cj > such that for every e > 
there is a N £ ¥1 such that 

F(|(i;,,T„„{.).....r,{.)t;,)| >e(^^^)") > 1 - e"'^" 
holds for all n > N and all normalized vectors ui, G C^. 

This gives us a first handle on the exponential decay. Note however, that the above corollary 
depends on z, but we need uniform estimates for some open arc of the unit circle. But using lemma 
14. m we can derive at least an initial scale estimate, valid for some fixed uq. 

Proposition 4.14. Let 9q n £¥l and suppose the group G^i^^ is non-compact and possesses no 
reducible subgroup of finite index. Then there exist o"0)7o S 1R+, a natural number uq > n and an 
open arc of the unit circle Is{Oq) centered around 6q with arc length 26 > 0, such that 

F (/or allOGls : \ {v, ,Tnje) T,{e) v,)\> e^""") > 1 - e""""" . 

Furthermore, the constants (7o,7o G F+ are independent of the normalized vectors vi, V2. 

Proof. According to corollarv l4.13| there exist 7, o" > such that for e > there is G N such that 
for m> N 

F (1 {v,,Tr,m ..... T,ie) v,)\ > e(^-^)"^) > 1 - e"-'" . 

If n > A^ we set no = n, otherwise we set hq = N. We can assume ^{Uno) = 1, see ()12p and lemma 
14.81 Hence, for e' > there exists k > such that 

F (3 i G {1, . . . , no} : \ai\ < k) < e' . 

By coroUarv 14.121 we have for arbitrary G T 



F ( I (^;, ,T,JB) .....TAe)v,)\> e^^"^^"" - r? ) > 1 - e--^"" - ^^\e - e.\ - e' , 



with a coefficient C(no) possibly depending on e' . First, we choose e > such that 7 — 2e > and 
set 7o = 7 — 2e. Then we choose r/ such that e^'>'~^-'"'° — r/ > e'^°'^°. The constants (Tq and e' are 
chosen in a way that assures e~°'°^° > e""""" + e' . Finally, a bound 5 > for the distance \6o — 9\ 
can be determined from the constant C(no) depending on k. □ 
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4.5. Thouless Formula. A Thouless formula relates the density of states with the Lyapunov ex- 
ponent. So we begin this section by defining the density of states for a disordered quantum walk in 
the usual way as a limit of measures corresponding to the finite dimensional approximations of the 
walk operator as defined in section [¥?2] 



Lemma 4.15. For every hounded and continuous function / : T — )• C ife have with probability one 
that 

(16) hm J ti{xNf{W^)) = ( {6o e, , f{W^)6o ^e,)+ {6^^e,, f{W^)5o )) ; 

where XN is defined in terms of projections onto the localized states Psi^ca '^s 

N N 

XN = ^ ^ ^ ^ ^Si(S)ea "I" -f'5_(jv+i)(8'e2 "I" ^5jv+i®ei =• ^ '] ^ ^ PSi(g)ea ~^ Pboundary 
l=-N a l=-~N a 

Remark 4.16. The reason why we do not restrict Wi^ only to the lattice sites —N to N, but include 
the states (5„(7v+i) ^2 o^nd (Jat+i ei in the definition of XN is that we want to compare it to the 
restricted operator Wi^{N), which is only unitary if we include these matrix elements (see eq. Q 



Proof. For every fixed / we can evaluate the trace on the left hand side of (jl6p in the standard basis 
of localized states 6i ej and get 

^ tT{xNf{W^)) =771^7T-rT(y^ i^l ® , f{W^)6l e, ) + tv{Pboundaryf{W^))) ■ 



4{N + 1) 4(iV + 1) ^ ' ^ ^ 

Since / is bounded by assumption, the second term on the right hand side vanishes when we make 
large. As in the self adjoint case one can now show that Birkhoff 's theorem applies to the random 
variables {Xi := {6i ej ,f(Wi^)6i (8) Sj)}, so for fixed / G C(T) we have a subset ilj C of 
measure one for which equation (jl6p holds. If we in particular look at a dense countable subset 
Co C C(T) we can consider the intersection of all the Qj with f £ Cq. Since this is a countable 
intersection of sets of full measure it has measure one as well and since it is dense in C(T), this 
finishes the proof. □ 

Equipped with this limit we can now define the density of states for a disordered quantum walk 
via the Riesz-Markov representation theorem: 

Definition 4.17 (Density of states). The density of states of a disordered quantum walk is the 
measure on the unit circle T defined by 

f(e) ^{dO) = ( (5o ® e, , f{W^)6o ® + (^0 ® , f{W^)6^ ® e, )) 

The unitary restriction W^{N) also defines a measure 'd^ on T via the Riesz-Markov representation 
theorem 

(17) ^ f{e) ^N{d9) = 4^^^^ tvifiWUN))) 

Later in this sections we connect these measures to the Lyapunov exponent in order to proof the 
Holder continuity of the integrated density of states as it is done in |Joy04| . Therefore we need to 
establish a relation between the unitary restriction Wuj{N) and the non-unitary restriction XnW^jXn- 
Such a relation is given by the following lemma, which is inspired by | Joy04| . 

Lemma 4.18. With the definitions from equation ^ we have that for every bounded and continuous 
function / : T — )• C 

J™ JtWTu (tr(/(W^.(iV))) - tr{xNf{W^)XN) = 

N-^oo 4(A -|- 1) 

holds. 
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Proof. Clearly it is enough to prove the above statement for functions with supremum norm equal to 
one. Since f(z) is continuous, for every e we can find a trigonometric polynomial Ps{z) of finite degree 
T(e), such that Ps{z) approximates f{z) uniformly on T up to e. Using the identity {Wuj{N))'^ = 
XAr(W^)"XiV for n G N and adding and subtracting p£{z) in the trace difference we get 

\tr{f{WUN))-XNf{W^)XN)\ = 

(18) = I tr{XN{Pe{W^) - Pe{W^)XN)) + tr((/ - Pe){WUN)) - XNif - Pe){W^)XN)\ 

< I tv{xN{Pe(.Wu.) - Pe(.W^)XN)\ + 8{N + l)e 

Where in the second step we used that Peiz) approximates f{z) together with the restriction by XN 
and thereby reduced the problem to the difference of powers of Wi^ and W^^. By induction we get 
the following identity for such differences 

k-l 

rk jirk\.. TTnrur tt-> \Tjy ^ 



XN 



i=0 



which holds also for negative powers k since we can then substitute the inverses by adjoints. Taking 
the trace we find for the summands on the right hand side 



< 



w,., - w,. 



In the last step we used that since W^^jXn is bounded and that since W and W differ at most in the 
two coins [/tv+i and C^-{Ar+i) their difference is of finite rank we can apply 



l^^lli < 



1 ll-^llop 



for A a compact and B a bounded operator [Sim05| . Substituting everything into equation (jlSp we 
get therefore 



1 



4(A^ + 1) 



\tr{f{WUN))-XNfiW^)XN)\ < 



2r(e) (W^-W^ 



4(Af + 1) 



where r(e) is the degree of Pe{z). Since — Wi^ is of finite rank its trace norm is bounded uniformly 
in N, the assertion is proven. □ 

The following proposition describes conditions on the eigenvalues of the unitary restriction Wi^{N). 

Proposition 4.19. For z & G we define the vectors 



(19) 



-1 



1 



ze 



-IT)'- 



and the spectral polynomial 

(20) pn{z) 



{(t>i{z),TN-...-T.NH{z)). 



Then z is an eigenvalue of the unitary restriction W^{N) of the Walk operator if and only 
if p]\f{z) = 0, hence all zeros of pi\f lie in T. 

Proof. Because of the structure of the restriction Wi^{N) the eigenvalue equation (Wi^{N) — z)(j) = Q 
implies, as in the infinite case, the following relations via transfer matrices between the components 
of the eigenvector (f) 

'cP{2l + 1) 
^<j){2l + 2), 

In addition we get the following relations from the boundary conditions imposed by the flip operations 



(21) 



Tdz) 



f(t){2i - 1; 



-2N - I) 



-2N) 



and 



z(P{2{N + I)) = e'" (^(2iV + l) , 
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which can be seen from ([TJ. Therefore the first and the last two components of (p have to be 
proportional to the normalized vectors 



<PUz) = ^I^J,^.) and ^^(^) = 7f( 



ze 



respectively. Because of relation (|2ip we can express 4>r{z) through <j)L{z) via a product of 2N + 1 
transfer matrices. Noting that the vector (pj^iz) as defined in equation ()19p is orthogonal to (pB^z) 
the vector 4>l{z) has to fulfill the desired relation 

{(t)i{z),TN ■ ...-T^NCpdz)) = Q , 
in order for z to be an eigenvalue. □ 

The Thouless formula now follows from (j20p . We utilize the following linear factor decomposition 
of the spectral polynomial 

4(Ar+l) 

(22) pr,{z)=z''''+^{<pi{z),TN-...-T^NMz))=CN- n (^-^''0, 

1=1 

where the 6i are the eigenvalues of the finite dimensional operator Wtjj{N) as defined in Section [4.21 
The coefficient Cat is determined from the transfer matrices T_7v up to T/y, more precisely, it is given 
by ^ 

Cm = e-^(^^+''") [ n a\ 

\l=-N / 

which follows from combinatorial considerations. Now, for \z\ ^ 1 we take the logarithm of the 
absolute value of (|22p and divide by 4(A^ + 1) to get 

2N + 1 \og{\{cPi{z),TM-...-T_Mz))\) ^ ^ log(|«d) , '^^'^ log(k-e^^'l) 

4(iV + l) 4(iV + l) ^2^^ 4(Ar + 1) + A{N + 1) ' 

Note that according to the arguments in section we may assume that all T; are well-defined, i.e. 
ai ^ for all I. Hence, by the law of large numbers the first term converges to the expectation of 
log(|a|) and with \z\ ^ 1 the second term can be written as 

^ 4iN + 1) = m + 1) " = J log{\z-e'^\)^NidX), 

with the density of states iD^ as defined in ()17p . The right hand side of this equation converges to 
the corresponding expression with density of states iD. In summary, we have the following expression 



log \z\ + — ^ E (log |a|) + / log(|z - e'^mdX) 



4(iV + l) ' 4(iV + l) 7V^oo2 

According to appendix [A] the second term on the left hand side of this equation converges almost 
surely to the Lyapunov exponent 7(2), which finally proves the Thouless formula 

(23) 7(^)= J log(|z-e^^|)^?(dA) + ^E(log|a|)-^log|z|, 

for 2; S C, l^l 7^ 1. As shown in appendix [Aj the Lyapunov exponent equals the limit of a deceasing 
family of subharmonic functions, since it is defined as the logarithm of the norm of some holomorphic 
matrix valued function. The right hand side of the above equation constitutes also a subharmonic 
function, and it thus follows from arguments which can be found in the classical work of Craig and 
Simon [CS83j . that the Thouless formula actually remains valid for z ^ T, i.e. \z\ = 1. 

The arguments used by Craig and Simon |CS83| also prove that the integrated density of states, 
as defined in the next section, is continuous on T as long as the Lyapunov exponent is non-negative 
on T. In fact, the integrated density of states is Holder continuous on T, which is proven in the next 
section. 
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4.6. Holder continuity of the Integrated Density of States. An important property we are 
going to use later on is the Holder continuity of the integrated density of states A/" on T, which we 
define for z € T in the following way 

/■arg(2) 

(24) Af{z) := / ^dX) 







This helps us to prove an estimate for the probability of cases where there is a value of the spectral 
parameter z £ T such that the denominator of the resolvent formula (|15p is small, called Wegner 
estimate, see section BTl The integrated density of states M inherits the Holder-continuity from the 
Lyapunov-exponent 'y{z), therefore we need the following lemma. 



Lemma 4.20. Let W^j he a disordered quantum walk satisfying the assumptions of theorem \3.1{ 
then the corresponding Lyapunov exponent j{z) is Holder continuous on T, i.e. there exists a finite 
constant C and a strictly positive number a such that 

\-f{z)--f{z')\<C-\z- z'\" ,z,z'€T. 

A proof of this statement for hamiltonian systems can be found in |CL90| . Lemma I4.2UI follows 
from a careful adaption of this proof to the setting considered in this paper. Although this is 
straightforward and contains no new ideas we give a self-contained proof of lemma [4.201 in appendix 

The Holder continuity of 7(z) in turn implies the same property for the integrated density of states 
A/", this is due to the fact that 7(2) and the density of states are related via the Thouless formula 



Proposition 4.21. Let W^j be a disordered quantum walk satisfying the assumptions of theorem \3.1[ 
then the corresponding integrated density of states J\f is Holder continuous on T, i.e. there exists a 
finite constant K and a strictly positive number /3 such that 

\Af{z) - M{z')\ < K\z - z']!^ ,z,z'eT. 

Proof. For a £ T and (j) £ [0,7r] we define the function Afa,<p{z) '■= ■^{z)x{e&<C:\avg{e-a.)\<<i>}iz) ^i^d 
apply the Cauchy transform for values of 2: G B satisfying | aig{z ■ d)\ < (j)/2. Note that according 
to |CS83| the integrated density of states M is continuous almost everywhere on T, so for |z| < 1 we 
can apply integration by parts to obtain 



T 



(25) = AA(ae^^)i^^ii^^ - AA(ae-^)i^^il^^^ + 

27ri 27ri 

+ j log(l - e'^z)^{dX) - j log(l - e'^z)^dX) . 

T |arg(6l-a)|>0 

In this formula we choose the logarithm as log(re*^) = logr -|- id with 9 S [— vr, vr) if | arga| < (f> 
and with 6 S [0, 27r) if | arga| > (j). Since M is almost everywhere continuous on T its Fourier series 
converges point wise almost everywhere to the exact value. The Cauchy transform projects the 
Fourier series of M to the positive powers, see e.g. |CMR06] . and since A/" is a real- valued function 
this implies the relation N'{z) = KN{z) + KJ\f{z) — KJ\f(0). Hence, the integrated density of states 
is Holder-continuous if the real part of each of the four terms in (|25p is Holder-continuous in the 
limit l^l — )• 1. This is obvious for the first two terms and the last term. The third term may be 
rewritten as 

log(l - e'^z)^dX) = j \og{\z - e'^\)^{dX) + i j arg(l - e'^z)^{dX) - iXo ■ 

T T 
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The constant Aq depends on the branch of the logarithm we have chosen and is real, hence, the real 
part of this expression is just the first integral, which is Holder-continuous by the Thouless formula 
and Lemma [4.201 □ 



4.7. Wegner Estimate. We proceed by proving a Wegner-type bound for one-dimensional dis- 
ordered quantum walks. We closely follow the arguments in |CKM87] and |DSS02| . In fact, all 
arguments more or less take over and there are only minor modifications to be made for the unitary 
case considered here. Let us state the result first and then adapt the proof. Recall the definition of 
leiOo) as the open arc of the unit circle centered around Oq with elements whose phases differ from 
the phase of only up to e. 

Proposition 4.22 (Wegner bound). Suppose that the single site coin distribution /i fulfills the con- 
dition stated in \3.1[ Then there exist for every < /3 < 1 and every a > a natural number 
uq = no{f3, a) and a = a{(3, o") > such that 

holds for all 9q and G N with N > no, and normalized vectors cf)^'^ G P2. 

The results of the last section are crucial for the proof of our Wegner-type bound. We need the 
following lemma, which is basically a consequence of the Holder continuity of the integrated density 
of states. It bounds the probability to find an eigenvalue of a localized eigenfunction in small sections 
of the unit circle, implying that the eigenvalues cannot clump together arbitrarily. Let Wa;(0, A^) 
refer to the unitary restriction of constructed in the last section, acting on the 2A^ + 3 lattice sites 
centered around the origin (which we indicate by 0). Recall that its image however is of dimension 
4(A^ + 1). 

Lemma 4.23. Suppose that the integrated density of states of the disordered walk operator is 
Holder continuous. Then there exist p > and C > such that for each Oq £T and for all < e < 1 
we have 

F 3 6* G leiOo) and if G C''^^+^) = 1 , such that 

{W^{0, N) -e)ip = and {\ip{-2N - 1)|^ + |(/j(2iV + 2)|2)l < e ) < A{N + l)Cef . 

Proof. Let E{0,N,9Q,e) denote the event whose probability, denoted by p, we would like to bound. 
Recall that it is defined with respect to the walk operator Wuj{0,N) acting on the 2N + 3 lattice 
sites centered around the origin. Now let E{k,N,9Q,e), A; G N, be the corresponding event for the 
unitary restriction Waj{k ■ 2{N + 1) , N) of W^j on the 2A^-|-3 lattice sites centered around k-2(N + 1), 
constructed in the same way like Wi^{0,N). Since the single site coin distribution is i.i.d on the 
different lattice sites, we get due to the ergodic properties of the system that 

p = lim / #{k G (-L,...,L) : iV, 6*0, e) occurs }. 

Now choose L G N and let ki, . . . , kj G {—L, . . . , +L} be such that the event E{ki,N, 9q, e) occurs, 
and let G T and ipi^.^ be the corresponding eigenvalue resp. eigenfunction of Wuj{ki ■ 2{N + 1), N). 
Let us extend these functions to the whole of (C('^^~^^y^(^+^) by setting them equal to zero outside 
their domain of definition, and denote these extensions by ip^^ . By construction, these functions form 
an orthonormal system. Consider now the unitary restriction VFtj(0, L ■ 2{N + 1) + A^), having image 
of dimension (2L + 1) • 4(A^ -|- 1). By the definition of the events E{ki, N, Oq^e) we have that 

11(1^^(0, L-2(iV + l))-^fc,)(^fcJ| < e yki. 

According to |ST85| . which is easily adapted to the unitary case, we thus have that the number 
of eigenvalues of Wt^{Q,L ■ 2{N + 1)) in the open arc Ie{Oo) is bounded from below by j. But this 
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number can be bounded using the integrated density of states N . We conclude that 
p = hm ^ #{k E (-L, ...,L) : iV, 6*0, e) occurs } 

< hm #{ eigenvalues of W^{{), L ■ 2{N + 1) + iV) in Ie{do) } 

= 4(iV + 1) ^lim^ (4(Ar + i|)(2L + l) *^ eigenvalues of W^O, L ■ 2{N + 1) + N) in /,(0o) } 

< 4(iV + l)AA(4(^o)) < Ci{N + l)eP. 

□ 

We proceed with the proof of the Wegner estimate. Again, we closely follow |CKM87] and [DSS02] 
and just briefly sketch the arguments. 



Proof of proposition \4-22 To begin with, note that it follows from corollarv 14.131 using routine ar- 



guments, that there exist a>0,(5>0, noGN such that we have for all G T, n > ?io, normalized 
vectors v, 

(26) IS (\\Tn{e) ■ . . . ■ TAe)v\\~^^ < e""!'^ . 

For a positive number t let = [tN^] + 1, where [.] indicates the integer part. Furthermore, let k 
be equal to Tai/26 and define e = exp(-cr(27V + 1)^). Following |CKM87| and |DSS02| , we see that 
after introducing the events 

A^{9,N) = [\\T_N+n^ie) • . . . • T_N{9)<P^\\ > e"^'} 
S,(^,iV) = {||r^i„^W.....r^i(e)</.^|| >e«^'} , 

we can upper bound the probability 

¥{3 6 e le{0o) ■■ \{4>^, TNie) • . . . • T_Nie)4>^ )\ < e) < (i) + {ii) + (m) + (iv) 

by four terms (i)-(if), which are defined as follows. 



(i) = F I 39 € IM : \{cl)^,TNie)-...-T^N(9)4>^)\<en f| {A^{e, N) n B^{9, N)) 



(ii) = F \A,{9^,N)nB^{9^,N)n \J (^A^{9,N) 

66/3,(60) 



(iii) = F \A^{9^,N)nB^{9o,N)n |J i^B^{9,N) 
\ 66/3,(60) 

{iv) = F{{A^{9„N)r) + F{{B^{9„N)r) 
Using Markov's inequality and equation (f26|) . we directly get 

{iv) < 2e-^^W' . 

Now assume that the event whose probability is (i) occurs. The transfer matrices Tj(6») as well as the 
vectors (f>^'^ correspond to a unitary restriction Wi^{N) of some walk operator. Recall that for any 
normal operator A, we have 

\\{A-z)-'\\ < — ^— 

infA6a(A) |A - Z\ 
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which together with lemma and remark [4.101 impHes that there exists 9 G /3e(^o) ^ a{Wi^{N)). 
Let (f G be the corresponding normahzed eigenfunction of Wi^{N). Since we assumed that 

both events A!i(9,N) and Bk{9,N) occur, we have that 

V,(-2iV-l) < ^e-^-^" 
^(2N + 2)) < ^e^^'^^'- 

Thus, we have that 

(i) < F 3 6* G heiOo) and G C'^(^+^) = 1 , such that 

(l^^(0, N) -9)ip = and M-2N - + |(^(2iV + 2))|2)l < g-i'^^' ) 

< 4(iV + l)C(max(3e,e-^'^^''))'' . 
To bound (ii) and (iii), we note that for 9 £T 



1 /det(C/i)(^- 



(27) r.(.)-T,.') = 

If the events ^^(^O) ^) ^ (0, A^)^ take place we have, if is large enough, 



for all 1] < 1. By applying again Markov's inequality and proceeding as in the proof of lemma l4.11| 

we get 

/njv+i \'n 

(ii) < 2^e-^-^' Yl E(||r,w||)'=-iE(||r,w-r,(e')ll)E(||r,(e')|ir'^+^-M 



k=l 



< 2ie-'"'^^ max(E(||r,(e)f)"'^+^E(||ri(0')ir)''^+^) (n7v + l)(3e)^ 



where we used (|a| + < |a|'' + and ^3) together with E < E(||ri(e')D- Now, by 

assumption, there exists rj > such that E (||ri(6'')||'') is finite. Hence, by choosing r small enough 
and applying the same procedure to {in), we can find some a2 > such that 

{ii) + {Hi) < 2 6^"^^" , 
and the result follows. □ 

We note shortly that the Wegner estimate can be extended to the following "uniform" version 
with respect to the quasi-energy ^ G T. Instead of bounding the probability to find an eigenvalue in 
some small interval, we would like to bound the probability that two independent disordered walk 
operators have one equal eigenvalue. Thus, in view of lemma consider the quantity 

F.,^ ( 3 ^ G T : I (0^ , r^(e) • . . . • T.Me)(t>^ ) | < e-'^^^A^+i)' 

and I (0^ , fjv(e) • . . . • f_jv(e)0^ )| < e-'^^^iv+i)") , 

where the transfer matrices Tj and Tj and the boundary conditions (f>^'^ and i^^'-^ are chosen in- 
dependently from another, and hence can be thought of as corresponding to two independent walk 
operators Wca^ and W^2> both acting on a 4(A^ + l)-dimensional Hilbert space. But if the event 
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occurs for some ^ G T, then we know that there exists an eigenvalue of W;^-^ in the interval of length 
g-CT{2Af+i) around 6. Using the fact that Wfj-^^ has 4(A^ + 1) eigenvalues, we can apply the union 
bound and end up with 

and I (0^ ,f;vW • . . . • T_^W0i^ )| < e"'^(2A.+i)^^ 

(28) < 4(iV + l)F(30 G /,-.(.^+.)4^o) : | , T^W • . . . • r_^(e)<^L )| < e-'^(^^+^)') 

< 4(iV + 1) e""^" 

< 3(2iV + l)e-"'(2^+i)' 
for some q' > 0. 

4.8. Multiscale Analysis. We now carry out a multiscale-analysis similar to the ones in |GKOH 
IKle08| in order to extend the initial length scale estimate to products of arbitrary length. The 
Wegner estimate is used to control the bad instances. 

Lemma 4.24 (MSA Induction step). Let A : T x — )• 1R,+ be family of positive, real valued random 
variables, indexed by an element of the unit circle. Assume that we can verify the following two 
estimates. 

1 There exist for any no G N some positive numbers 70, cr and a subset /(no) of the unit circle 
such that 

F ( for all 9 G /(no) : X{9,oj) > e^°"°) > 1 - e"^"" . 

2 Given a natural number m and the product 

m 

1=1 

where the factors are understood as being independent and identically distributed (i.i.d), we 
have for all < 6 < 1 the estimate 

(29) F ( 3 G /(no) : A'"(^, w) < e"^"** and A'"(0,w) < e"^"**) < Sme"'^'"' , 
for two i.i.d. products X'^{6,io) and X'^{6,uj). 

k 

Then we have for all natural numbers k, 1 < a < 2, ^ + 1 < a, length scales = 

(30) F ( 3 G /(no) : A"'- (^, w) < e^^"* and X''^{e,ui) < 6^^="*) < e-'""^ 
holds with 7o > 7fc ^ • 

Proof. We prove the lemma by induction. So assuming that (|30p is fulfilled for some /c G N, we use 
assumption 2 to prove that is also holds for k + 1, with a slightly smaller 7^+1, which can nevertheless 
be lower bounded by 7o/2. We find that this requires that, depending on our choice of a and ^, the 
initial scale no has to be large enough. The result then follows by applying the first induction step 
to a suitable initial scale no, which can be constructed using assumption 1. 

We divide the proof of the induction step into two parts, the deterministic and the probabilistic 
estimate. In the deterministic part we make assumptions about the value of the product of the 
random variables A"''= that helps us to estimate a new growth rate 7^+1. 

In the probabilistic part we bound the probability that these assumptions are violated or in other 
words we bound the probability that the new growing rate is smaller then 7^+1 as estimated in the 
deterministic part. 

So let us begin by stating the two assumptions we use in our deterministic estimate. 

Al Only n^ of the n'^~^ factors in the product of the X"'''{9,oj) grow with a rate smaller then 
7/j. Where < /5 < 1 is chosen such that a > 1 + /3 
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A2 None of the bad instances from assumption Al is smaller than e~~^"'=. 
Now let us find a lower bound 7^+1 on the growth rate of X'^'^+i. Taking the logarithm we get 

(31) ln(X"'=+i {e,uD)=Y, HX""' (e, Ui)) > - ^hkUk - ^ninl 

i=l 

where we used the assumptions Al and A2 in the second step of the inequality. In order to lower 
bound the new growth rate we define 7^+1 by the right hand side of (|3ip divided by n^. This leads 
to the recursion relation 

(32) 7fc+i = {\-n^ )lk - yn;, > 7fc - ip^n^ , 

because b < \ and 7^ < 70 by assumption. It remains to show that still 70 > 7fc+i > "1^ holds. Let 
us start with the lower bound. Considering (|32p we find 

(33) 7fc+i - 70 = X^(7i+i - 7O > -^70 X] A''^''^ ■ 

1=0 1=0 

Using nfc = tiq'' , and setting q = nj^'^"" we have to show that 

k 

(34) ^ + ^^"'<i 

1=1 

to guarantee 7^+1 > Since the summands are positive and q < 1 holds, because 1 + /? < a, we 
can bound the finite sum by the infinite series and use 

< Inserting this into ([33]) . 

evaluating the geometric series and exploiting q < g"~^ we get the condition 

(35) < 1 

which can be achieved for uq large enough since 1 + /3 < a. By assumption this proves 7/0+1 > 
So we are left to prove that 7^ is a decaying sequence. From 

,l+/3-a„,, , 70^^<5+/3-cv>, 

2 " 

we conclude 7^ > 7^4.1 and therefore 70 > 7fc+i by assumption. This finishes the deterministic 
estimate. 

Let us now proceed with the probabilistic part, i.e. we have to bound the probability that either 
assumption Al or assumption A2 is violated for both products X"'=+i(^,a;) and X'^''+^ (9 , uj) . Let 
us denote the event that assumption Ai is violated for X^''+'^ {9 , uj) by CAi[X"''+i (^, ;^)]. We can 
estimate 

F{39 e I{no) : Assumptions Al or A2 are violated for (6^, w) and (6^, w)) 

= 2F(30G/(no) : CA1[X"*+- (0, w)] and CA2[X"'=+^ (0, w)]) 
+ F(30G/(no) : CA1[X"'=+-(0,w)] and CA1[X"'=+^(0, (!>)]) 
+ F{39e /(no) : CA2[X"^+- (0, w)] and CA2[X"'=+^ (0, w)]) 
< 3F {3 9 eI{no) : CA1[X"'=+- (0, w)]) 
+ F{39e I{n„) : CA2[X"'=+- (0, w)] and CA2[X"^+^ (0, (!>)]) 



(36) 7fc+i -7fc = -(n^ ^ 7fc + v'^fc )<0 



Now assume Al is violated, i.e. more than factors of the product X'^{9^ijj) grow with a rate 

smaller then 7^, or put differently, more than — pairs of products grow with a rate smaller then 
7fc. Hence we can use the induction hypotheses (|29p to bound 

3F(3eG/(no) : CA1[X"''+^ (0, lj)]) 

< 3F(3 6'G/(no) : CA1[A"'= (0, a;)] and CA1[X"'= (0, lD)]) ^ < Se"'""^-^' 



fc 2 
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From the requirement (j29p we get that 

F(3^G/(no) : CA2[X"'--+^ (0, w)] and CA2[X"'=+- (0, w)]) < Snfce"'""^. 
Thus we can combine these two estimates in order to get 

¥{39 £ /(no) : Assumptions Al or A2 are violated for X"-''+^{e,uj) and X^^^+i (6', w)) 

provided that ^ + /3 > a • ^ and 5 > a ■ and uq big enough. 

Together with the requirement f3 < a — 1 originating from the deterministic estimate, we see that 
the induction step can be carried out if,^</3<a— 1 and ^ < 6 /a. But ^ + 1 < a holds by 
assumption, and thus we can find /3 and 6 accordingly. □ 

4.9. Dynamical Localization. 

Lemma 4.25. Assuming the conditions of theorem \3. 1\ are fulfilled, then there exist positive constants 
Ci; (^2; < ^ < 1 such that we have 



where is the spectral measure of the walk operator corresponding to sites x, y £lL. 

Remark 4.26. The right hand side of this inequality clearly constitutes a function 

which decays faster than any polynomial in n, and together with the arguments given at the beginning 
of section \4-^ the preceding lemma thus also proves theorem \3.1\ 

Proof. We show that if the conditions of theorem 13.11 are fulfilled, then there exists for every ^ G T 
an open arc I5, 5 > 0, around 9 and appropriate constants Ci, C2, < ^ < 1 such that 



for large enough ring sizes L. The result then follows by compactness of T. 

For the rest of this section, we assume without loss of generality that x > y. Using lemma and 
the theory of Cauchy transforms (see section 2]2|) , we can apply the upper bound 



els : 2\{<^+,T^_i{z)...Ty{z)^_)\ < - 

K 



Hence, we have to lower bound the growth rate of | ,Tx-iiz) . . .Ty(z)^-)\. This is done by 
applying the MSA lemma. 

We first choose anl<a<2,0<^<l, and an initial length scale hq G N respecting the 
requirements of lemma 14.241 such that there exists A; E N with the property that the /c-induction 
step applied to uq provides us with a length scale i = Uq fulfilling 2i + 2 < \x — y\ < , where 
is identified with is integer component. Note that this is always possible by suitably choosing no, a 
and ^, if |x — y| is big enough (the other cases are of course not important if we look at the decay 
rate, and are covered by a constant). Then, we can factor out from | , T^^iiz) . . . Ty{z)^^ )| a box 
of length i, either starting at x and stretching out to the right or starting at y and stretching out to 
the left, i.e. 



{<!>+ , r,_i(.) . . . Ty{z)<i>. )! = !($+, T,_i(.) . . . r,_i_K.)|._ )||(|.+ , r,_2-K.) . . . Ty{z)<^_ )\ , 

I (^'+ ,T,_i(.) . . . Ty(z)<!>^ )| = I ($+ , . . . )||(|.+ , Ty+eiz) ■ ■ ■ Tj,{.)$_ )| . 
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Both of these boxes can again be factored into i.i.d. boxes of length no, i.e. 

I , r,_i{.) . . . r,_i_^(.)l>_ )| = J] I (1.^ , T„„{.) . . . T^{z)¥_ ) I 

with some collection of normalized vectors {$!|-} in the plane Pi, and the same can be achieved 
for I ($_|_ , ra;_i(2) . . . Tj^_l_f+i(2;)$_ )|. Now by assumption, the conditions of corollary 14.131 hold for 
the group of transfer matrices associated to any G T. Hence we may apply proposition 14.141 to 
construct an open arc J^, 5 > and get an initial length scale estimate valid for the length no with 
constants 70 (for the scalar product) and u/ (for the probability). Thus, assumption 1 of lemma 11.241 
is verified for the random variables 

|(^>V,r„oW---T^iW^-)l ■ 

The verification of assumption 2 of lemma 14.241 is done using the Wegner estimate and was carried 
out in the discussion preceding equation (j28p . with decay constant tuy (bound on the probability). 
Note that we can choose a = min((Tiy, <7/), and both the Wegner estimate and the initial scale 
estimate are still valid for a as the decay constant in the bound of the probability. Hence, both 
assumptions needed for carrying out the MSA are verified. Next, define the event A{Is,tj <Z as 
the set of those lo for which either the box 

|($+,T,_i(^)...T,_i„K^)^->| 

or the box 

is bigger or equal to e^^, 7 = ^70 for all 9 € 1$. Here, 70 is the deterministic decay constant we get 
from the initial scale estimate. As assumptions 1 and 2 of lemma 14.241 are verified, we have that 

(37) F{A{Is,i))>l-e^'''' . 

Assuming that A{Ig,i) occurs, and reshuffling, as well as renaming indices, we end up with 

1 



(38) 
(39) 



"^'^rK/^) < lim TTK ■ X 



\Puj,n\ 



els : e^^2|($+,2|,„^|„i„^(.)...ri(.)$_)| < 



lim TTK ■ A 



els : 2|($+,2|,_,y|„i_,{.)...ri{.)$_)| < - 



is, by lemma [4^9 



with ^± being vectors in the plane P2. But |^2 ,T\x-y\-i-£i^) ■ ■ - Tiiz)^. 
and remark r4.10| the resolvent of some unitary walk operator restricted to a box of length |x — ?/| 
It then follows from equation ^ that 

1 



(40) hm TTK ■ A 



els : 2|($+,2|,_,y|_i_^(.)...ri(.)^>_)| < 



Now, if we combine equations ([37]), ()38p and ()40p we can estimate 



E 



A{Is/) 



Fide) ip^'im + 



ZAiIs,£) 



< 1 . 



< e 



+ e~ 



under the assumption on |a; — y|. But since we have chosen £ and a such that |a; — y| < it follows 
that 



(41) 



^{\P:ViIs)) <C{e- 
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where the constant C > is chosen to cope with the cases where |x — y| is not "big enough". Given 
the requirements of theorem I3.H the derivations in the preceding sections assure us that we can 
always find 7, cr > and ^, a such that equation ()4ip is true with < ^ < 1, and the proof is 
complete. □ 

What remains to be proven is the decay of the eigenfunctions. We first prove the statement for 
the eigenfunctions of the finite restrictions W^{N). This can be derived from our bounds on the 
time evolution of initially localized particles. We use a variant of Wiener's theorem for measures on 
the unit circle, proven in the appendix, see theorem ID. II Using in addition that the matrix elements 
of some eigenprojector P'^q^ of W^{N) at the eigenvalue 6 are given by the value of the spectral 
measure at the point ^, 



{5y ® e. 



we arrive at 



6lG(7(W„(Af)) 



lim „ 

r^oo r + 



1 ^ 

—y 



i=0 
T 



< lim V I ei , W^(iV)*4 ® ) p 

t=0 



< sup I {6y (E>ei,Wu; {NYS^ ) | 

ieIN 

where we used the Cauchy-Schwarz inequality in the last step and a{W;^{N)) denotes the spectrum 
of Wi^{N). Now, we have by lemma [4.251 that we can bound the decay of the expectation over all 
configurations of the spectral measure by £, which implies 



El E E ^(1^ 



y\) ' \ {6y ^ ei , P^QySa: 



< 00 



since C decays faster than any polynomial. This in turn implies for every 9 G cr{Wi^{N)) and any 
w G the existence of a random variable with finite expectation Ag{uj), such that 



{6y (g) e 



« ' ^{9} ^x 



)|2 < Ae{oj)C{\x-y\), 



ensuring the decay of eigenfunctions of Wuj{N) with probability one. To apply the argument for 
the eigenfunctions of W^, defined on the infinite lattice, we have to show that the spectral measures 
Pu% converge to pZ'^ of W^, if — )• 00. But this follows from equation ([3]) and the fact that the 
polynomials are dense in the space of continuous functions on the unit circle. 



Appendix A. Products of random matrices 

In this section we gather some parts of random matrix theory. Our aim is to analyze the behavior 
of random products of matrices, in particular we are interested in the growth of quantities like 

\\gn ■ . . . ■ gi ■ v\\ and \ {w , g-n ■ ■ ■ ■ ■ gi ■ v)\ 

where v and w are normalized vectors and the g^ are complex valued matrices drawn independent 
and identically distributed according to some measure n on GL(K,/c) where K = IR or C In fact, 
the measure ji depends on a parameter 2; G T, which we take into account by writing n^, hence, the 
matrices gi also depend on z but to simplify notation we omit this dependence. 

If K = IR all results derived in this appendix are well-known, see e.g. |BL85| ICL90| . however, for 
K = C the literature is less exhaustive. So, in this appendix we adapt these results to the case 
K = C Many statements can be translated directly from IR, to C by observing C = IR^. This 
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approach is similar to the analysis carried out in |BHJ03| . More precisely, we identify C with 
via the isomorphism 

,2 ■ ' ■ - 1^ 



(42) k : <C ^ ^\ h{x + iy) = \^^y x,yGlR. 

This mapping naturally extends to and induces a homomorphism ^ mapping linear operators 
on to linear operators on via the relation 

^ : Mat(C, k) Mat{K, 2k) , ^'(M) = A o M o . 

It is easy to see that both A and ^ are R-linear and since A is an isomorphism ^ is injective, however, 
it is not surjective. Moreover, ^' is a homomorphism with respect to matrix multiplication, hence a 
subgroup G of GL(C, k) is mapped isomorphically to a subgroup ^{G) of GL(1R,, 2k). The mappings 
A and ^ have the following useful properties: 

Proposition A.l. The mappings A and ^ preserve the euclidean norm, i.e. 

||M-^;||^ = ll^(M) • A(t;)||j^, VM G Mat{(C,k),v G , 

where the norm is given by the euclidean norm in C'^' respectively R^'^. Moreover, we have that the 
map ^ preserves the modulus of the determinant, i.e. 

|det(M)| = |det(tp'(M))| , VM G Mat(C, k) . 

Proof. The first statement of the lemma follows easily from the fact that 

\\^{M).A{v)\\^ = \\A{M.v)\\^ 

for ah M G Mat(C, k) and v G C'^. 

Since the modulus of the determinant of a matrix M is given by the product of its singular values, 
it it sufficient to prove that ^ preserves the singular values of M. This follows from the equations 

*(M • M*) = *(M) • ^'(M)* , VM G Mat(C, k) 

and 

A(M • v) = ^'(M) • A(v) , VM G Mat(C, k) , v £ & 
together with the invertibility of A. □ 

The transfer matrices T we consider satisfy |det(T)| = 1, therefore we are interested in ^(T) 
where T is an element of the group of matrices with determinant of modulus one, which is denoted 

by 

SLT(fc) = T X SL(C, k) = {M G Mat(C, k) : |det(M)| = 1} . 

Remark A. 2. Since the map ^ is continuous we have that det(.) o is a continuous map from 
SLt(A:) to { — 1,1}. The group SL(C A;) as well as SLt(^) is simply connected. Therefore, the 
function det(.) o ^ is constant on SL^^^lk) and since det{\I^{l(^k)) = det(lj^2fe) = 1 it is equal to 1. 

Our aim is to prove that, under certain circumstances, the Lyapunov exponent, as defined sub- 
sequently, is positive and Holder-continuous on T. To this end it is helpful to consider the group 
generated by the support of the measure /i^, which is defined in the following way: 

Definition A. 3. Let {fiz}z£T a family of measures on GL(K,A;), the general linear group on the 
'K-vector space of dimension k. 

i) The Lyapunov exponent associated with /i^ is the limit 

7(z) := lim -E(||5n-...-9i||) , 

n^oo n 

where the expectation value is with respect to the n-fold convolution /i" of the measure jiz. 

ii) The support of the measure is defined as the set 

supp(/x^) := {M G GL(K,A;) : iiz{Be{M)) > OVe > 0} , 
with Bs{M) denoting the sphere of operator norm radius e around M G GL(K, A;). 
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iii) The smallest closed subgroup o/GL(K, A;) containing the elements of supp{fiz) is denoted by 

il^z)- 

To such a measure we may associate a map acting on functions / from the projective space 
PK*^ = {v £ K'^ : ||f|| = 1} to the real numbers in the fohowing way 



(43) {Rzf){x) := J f{gx)^z{dg) for bounded measurable / : P'M} 

From now on we adopt the notation x for elements in PK'^, whereas x denotes the elements of K'^. 
In the following, we restrict the definition of Rz to functions / from a set defined with the help 
of a distance 8 on P'K.^ and a > in the following sense 

fcr ^ l/(^)-/(^)l . 

f mo^U) := sup — — < oo . 

x^y <5(x,y)" 

If we define a norm on via ||/|| = ||/|| + niaif) it becomes a Banach space (note that / 
is assumed to be bounded), and denoting the set of bounded operators on Ca by B{Ca) we have 

This leads us to the notion of an invariant measure Vz corresponding to /i^. 

Definition A. 4. Let {liz^zel a family of measures on GL(K, A;). An invariant measure Vz for 
[iz is a measure fulfilling 

f{gx)^^{dg)vz{dx) = / f{x)u^{dx) 



for any bounded and measurable function f : P'K'' — )• R. The operator Nz G B{Ca) defined via 

{Nz{f)){x) := J f{y)uz{dy) 
is called the projection onto constant functions. 

Remark A. 5. According to Prokhorov's Theorem the space of probability measures MiPR^) consti- 
tutes a sequentially compact set, i.e. each sequence contains a weakly-converging subsequence. This 
fact can be used to give an expression for an invariant measure for a measure fx. Let 6 be an arbitrary 
probability measure on , then 

^ n— 1 

n ^-^ 

k=0 

defines a sequence of probability measures on R^' and it is easy to check that the limit of any weakly- 
convergent subsequence is an invariant measure for /i. 

While invariant measures Uz are hard to construct explicitly, their existence and certain properties 
can be established by analyzing {fJ-z)- In fact, many properties of 7(z) can be established through 
analysis of the corresponding invariant measure Vz- We start with the connection between the real 
invariant measure and its complex analogue. 

Remark A. 6. Let fi be a measure on GL(C, k) and denote the corresponding measure on GL(1R, 2k) 
obtained via the isomorphism A by An invariant measure for ;U]r induces an invariant measure 
V for /X via 

v{x) = U k{ex)) , V X c pc^ 

Hence, if there is an invariant measure for the real embedding of fi then there is an invariant measure 
for ji itself. 

In particular, for K = IR there exists a unique continuous invariant measure Vz if {l^z) is strongly 
irreducible and non-compact |BL85| IGL90| . 

Definition A. 7. Let G be a subgroup o/GL(K,A;) with /c G N. We call G strongly irreducible over 
K if there is no finite sequence of nontrivial subspaces Si, . . . ,Sn C K'^ such that G {Si U . . . U Sn) C 

SiU...USn. 
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Remark A. 8. A subgroup G C GL(K, k) is strongly irreducible over IK if and only if all subgroups 
of G with finite index are irreducible over K. 

The following two propositions translate non-compactness and strong irreducibility from R to C 

Proposition A. 9. A subgroup G of GL{(C,k) is compact if and only if^{G) is compact. 

Proof. The map ^ is continuous, hence ^{G) is compact if G is compact. On the other hand, 
we can invert \I' continuously on the set ^'(GL(C, /c)). As a consequence G is compact if ^(G) is 
compact. □ 

We need to identify circumstances under which the subgroup ^{G) is strongly irreducible. For G 
this means that all subgroups of finite index are irreducible over R, i.e. the only R-linear invariant 
subspaces of C'^ are trivial. The following proposition simplifies the analysis in the sense, that it is 
sufficient to check strong irreducibility for G over C in certain cases. We prove later (see the proof 
of lemma PV.14p that these are exactly the cases of interest to us. 

Proposition A. 10. Let G be a subgroup o/GL(C,/c) and denote the direct product of G and T by 
Gj = T X G. The group ^(Gj) is strongly irreducible over H if and only if Gj is strongly irreducible 
over C 

Proof. Clearly, if ^{Gf) is strongly irreducible over R, then G^ is strongly irreducible over C Now, 
suppose there is a finite sequence of nontrivial subspaces Si, . . . ,Sn C such that ^(Gt) Si U 
. . . U S'n C U . . . U 5„. The sets Pi = k-^{Si) C C'' satisfy Gt ^'i U . . . U C Pi U . . . U P„ and 
each Pi is an R-linear subspace of . But since Gf contains all operators e***^!, cj) G R each Pi is 
also a C-linear subspace of C'^, hence Gf is not strongly irreducible over C □ 

The next lemma summarizes the results we wish to translate from the case K = R to K = C For 
a proof of the first part we refer to |Fiir63| or |CL90| proposition IV. 4. 6, the second part can also be 
found in |CL90| propositions I.V.4.4 and IV. 4. 5, and the third part is a variation of theorem V.6.2. 
in |BL85| . First, we define the term C-integrability of a measure ^. 

Definition A. 11. Let fi be a measure on GL(K, A;) with k £ ¥1. We call fi (-integrable, if there 
exists C > such that 

(44) j \\g\\^ n{dg) < oo . 

Remark A. 12. For unitary matrices U and z G T the norm of the transfer matrix 
^^iU) = -(^ l), wUh U=(l 



a 



can by bounded by \\tz{U)\\ < 2/\a\. Therefore, IE ^||P||^j is finite for some ^ > i/ IE (|o| 
finite. 



IS 



Lemma A. 13. Let {fiz}z£T be a family of (-integrable measures on SL(R, /c). If the groups {^z) C 
SL(R, k) are non-compact and strongly irreducible over R, then 

i) The Lyapunov exponent satisfies 7(2) > for all z G T. For all vectors v ^ in R'^ 



lim -log\\gn ■ . . . ■ gi ■ v\\ = j{z) 

n— >oo n 



with probability one. 
ii) Defining the function 



$,(x) := [ log ^p^//, (dp) 
we have the relation 7(2:) = Nz{^z) and for fixed z £T the function ^z is continuous. 
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iii) There exists a > such that for any v ^ in M!' and e > there is N with 

F{\\gn...gM\ >e"(^-^)) >l-e-'^" 

for all n> N . 

Given a measure ^ on a subgroup G C GL(C,A;) we can extend this measure by an arbitrary 
measure u on T, i.e. we consider the product measure fiy = fiy. v onT x G. For arbitrary g £ G and 
gi(/> g TP it is clear that | |e*''^ x g\\ = \\g\\. Hence, we can simulate the random variable \\gn ■ ■ •^I'yll, 
where the gi are drawn according to fx, by drawing the gi according to iJ,iy. But this means we may 
equivalently consider the group (fj-iy) instead of {^) when analyzing the norm growth of the matrix 
products. If we choose the uniform distribution on T as we have (fi^) = T(/i) since supp(z^) = T 
in this case. 

The next lemma is the desired generalization of lA. 131 for complex valued matrices. 

Lemma A. 14. Let {fiz}zeT be a family of C,-integrahle measures on S\jj{k). If the group {^z) C 
SL']r(/i;) are non-compact and strongly irreducible over C, then 

i) The Lyapunov exponent satisfies 7(2) > for all z G T. For all vectors v in C^' 



lim -log\\gn ■ . . . ■ gi ■ v\\ = j{z) 



n— >oo n 



with probability one. 
ii) Defining the function 



:= /"log^^/i,(d5) 

J IfII 



we have the relation ^{z) = Nz{^z) and for fixed z £T the function is continuous. 
iii) There exists a > such that for any v ^ in and e > there is N £ ¥1 with 

F{\\gn...gM\ >e"(^-^)) >l-e-'^" 

for all n> N . 

Proof. For the most part, the lemma is just a translation of lemma lA.131 from C*^ to IR^*^ using the 
propositions of this section. It remains to prove that without loss of generality we may consider 
the group T(/i) instead of such that we can apply corollary lA.lOl But, as already mentioned, 
the norm of an element gj x gQ £ T x G only depends on the factor go- Moreover, and go are 
independently drawn according to the product measure ^xv where v is the uniform distribution on 
T. Therefore, the cumulative distribution function Fx{x) = F {X < x) for the two random variables 
X = \ \gn . . . giv\\ with gi either given by go or gj x gQ coincide and the lemma is proven. □ 

Equipped with these basic results we now consider the space and prove certain regularity 
properties of the invariant measure z^^ and later on the Holder-continuity of 'y{z). Still, these results 
are known for K = R, but their translation to the complex domain is less obvious. The first lemma 
we now prove gives the invariant measure i^^ an interpretation as the probability measure of a random 
variable on PC^. Let us denote the n-fold left product gi gn with gi £ supp(/i^) by m". 

Lemma A. 15. Let fi^ be a C,-integrable measure on GL(C2); with |det(5f)| = 1 for g G supp(/i2) 
and suppose that {fi^) is non-compact and strongly irreducible. It follows that m" • converges 
almost surely to a rank one matrix and its direction is a random variable which is distributed according 
to the invariant measure Vz corresponding to fiz- 

Proof. It is clear from lemma IA.14I (i) that if the Lyapunov exponent is positive then the biggest 
singular value of the product m" has to grow exponentially with probability one. Since | det(m")| = 
1, the second singular value is equal to the inverse of the first one, and hence decreases exponentially. 
Hence m" • ||m"||~^ converges almost surely to a rank one matrix. Let mz{oj) denote any such limit 
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point. Since (^u^) is strongly irreducible, the invariant measure is continuous and hence the 
equation 

j(x)v^{m:,(bj)dx) = f f ( TT-^T^^^] ^z[dx) 



\mz{'jj)x\ 

defines a measure mz{u))i^ on PC^. If m" • ||m"||~^z^ is defined accordingly, but with mz{oj) replaced 
with m" • ||m"||~"'^, then the measures m" • ||m"||~-^i^ converge weakly to mz{oj)v. Moreover, by 
lemma II. 2.1 in |BL85| . the measures m" • ||m"||~^z^ converge to some probability measure on PC^, 
which according to the above considerations is equal to a Dirac point measure Hence we have 

mz{oj)v = 5d{ui)- This proves the assertion. □ 

Lemma A. 16 (Iwasawa Decomposition). Let M G GL(C, k) with A; G N. Then there exist a unitary 
matrix U{M) and a lower triangular matrix s{M), with positive diagonal entries such that 

M = s{M) ■ U{M) 

Moreover, this decomposition is unique and we have that 

MM*ei _ s{M)ei 
||M*ei||2 ~ (s(M)ei ,ei) ' 
where ei, . . . , e„ denotes the standard basis of C*. 

Proof. The existence, as well as the uniqueness of the decomposition, follows by applying the Gram- 
Schmidt procedure to the vectors M*ei, . . . , M*en- The second statement is easily seen after noting 
that s(M)*ei = (s(M)ei ,ei )ei. □ 

Lemma A. 17. Let fi^ be a C,-integrahle measure on GL(C,2); with |det(5f)| = 1 for g G supp(/i2) 
and suppose that (fiz) is non-compact and strongly irreducible. Then there exists an a > such that 

\s{V*m^V)e^\ 



supE 
V 



< oo , 



where V runs over the unitary subgroup o/GL(C,2) and s{V*m^V) denotes the lower triangular 
matrix in the Iwasawa decomposition of V*m^V . 

Proof. To shorten notation within this proof let us define for a matrix m 

sirfi) 

(45) rim) := — r and rjim) = r(m)ei — ei 

(s(m)ei ,ei ) 

Therefore our goal is to bound the expectation value of ||T(y*m"y)ei||" from above for some a > 0. 

Let us now study the Iwasawa decomposition of a product of two matrices mi and m2 . Applying 
the decomposition either to the full product or subsequently we get 

mim2 = s{mim2)U{mim2) 

mim2 = s{mi)U{mi)m2 = s{mi)s{U{mi)m2)U{U{mi)m2) 
which implies by the uniqueness of the decomposition 

(46) s{mim2) = s[mi)s{U{mi)m2) 
Inserting this identity into our definition of r] we find 

n-l 

r/(y*m^y) = r]{V*miV) + ^T{V*m^^V)r]{U{V*m^V)V*Yk+iV) 

k=l 

for a left product m" of matrices mj. Returning to the quantity of interest we find 
\\T{V*m^V)ei\\ = 

n-l 

= ||r?(T/*m^y) + ei|| < ||T(y*miT/ei)|| +J2\\'^i^*'^^^)viU{V*m^V)V*mk+iV) 

k=l 
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The first term on the right hand side is bounded from above by one so we are left to study the 
remaining summands. First note that rj^m) is by construction orthogonal to ei , so in order to upper 
bound it, we have only to consider vectors along the 62 direction. Rearranging the the expression a 
little bit and abbreviating U (V* m}^V)mkj^i by x we get 



(47) 



s{V*m^V)ri{x)\ 



Ux)\\ 



Mx)\\ 



{s{V*m^V)ei ,ei 



\\s{V*m^V)e2\\ 
{s{V*m^V)ei ,ei)\ 



The expectation of the first factor can be bounded by properties of the Iwasawa Decomposition (see 
lemma P\..16p and the determinant condition by 



(48) 



E(||r?(x)|r)<E 



|x*ei||' 



<]S{\\x\\\\x*-Y)<]S{\\mk+,\n 



Due to the C-integrability the expectation value on the right hand side is finite for some a > and 
independent of k. We now consider the second factor of equation (j47p . Since ||s(x)e2|| = {s{x)e2 , 62 ) 
and the lower triangular form of s{x) for every x we find 



(49) 



E 



\s{V*m^V)e^ 



< E 



dets{V*m^V )Y 
\\V*m^VeA\ 



< E 



Since 7(2:) is positive there are constants C > and < /3 < 1 such that this expectation value can 
be bounded by (7/3*^. Putting everything together we get 



(50) supE 
V 



\\s{V*m^V)e^\\ 
(s(y*m^y)ei ,ei 



<E(||77(y*m,y)ir) + ^E 



E(||mJ 



which is finite for some < a < 1, because the infinite series is bounded from above by a converging 
geometric series. 

□ 

This leads us to the following corollary, which provides us with an estimate on the regularity 
properties of the invariant measure. 

Corollary A. 18. Let fi^ be a Q-integrahle measure on GL(C,2), with \ (lei{g)\ = 1 for g G supp(^2) 
and suppose that (^^) is non-compact and strongly irreducible. Denote its invariant measure by v^- 
Then there exists an a > such that 



sup 



{x,y)\ 



Vz{dx) < 00 . 



Proof. Let ai,2('T') be the singular values of the n-fold product m" and note that because we have 



det{g)\ = 1 for 5 € supp{fj,z), it follows that 02 (n) 



ai[n) 



Since (/U^) is non-compact and 



strongly irreducible, we have that ai{n) — t- 00 for n — t- 00, and moreover, the direction of the rank 



one matrix lim m"\\m^ 

n— >oo 



n||-l 



is distributed according to Uz- Let 



be the polar decomposition of m", with K'^ being unitary matrices. It follows from the above 
arguments that the distribution of the limit lim U^ei is again given by i^z- Now, since for any unit 

n— >oo 

vector u E C^, 

(m^(m^)*ei ,n) = aiinf {U^eud ) (d , C/>) + 0{a^\n)) + 0(1) 
we have by lemma [A. 161 that 

s(m")ei D 
hm —— = 

n^-oo (s(m^)ei ,ei) {D,ei) 
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where D is a random variable with distribution Vz- Next, let y G P<C^ be some unit vector and let 
V be defined hy y = V*ei. If we replace each matrix g in the n-fold product m" by VgV* , then all 
arguments remain valid up to replacing m" by Vm^V* and D by VD. It follows that 

s(VmW)ei VD 
lim — 



n- 



>c» {s{Vm^V)ei,ei) {D ,y] 



Since the distribution of D is u^, we have that 

J \\{x,y)\J ' - \[{s{V*m^V)e,,e,)\ J 

which is finite for some a > by lemma P\.. 171 □ 

Finally, we give a proof of corollarv 14.131 To be more precise, we prove the following more general 
lemma, which implies [4.131 when taking into account the properties of the transfer matrices. 

Corollary A. 19. Let ^ he a C,-integrahle measure on GL(C, 2), with \ det{g)\ = 1 forg S supp(^) and 
suppose that (/x^) is non-compact, strongly irreducible and contractive. Then there exists 7, a, Eq > 
such that for every Eq > ^ > there is a N € ¥1 such that 

F(|(t;,,r„„{.).....r,{.)t;,)| >e(^-^)") > 1 - e"'^" 

holds for all n > N and all normalized vectors vi, V2 £ C'^ ■ 

For the proof of this corollary we follow along the lines of [DSS02| . The actual proof for corollary 
lA. 191 relies on the following result. 

Lemma A. 20. Let ^ he a Q-integrahle measure on GL(C,2), with \ (iet{g)\ = 1 for g S supp(^) and 
suppose that (nz) is non-compact, strongly irreducible and contractive. Then there is a Eq > such 
that for all < e < Eq there is a 6 > and E N such that for all n > N and y G P<C^ 

(51) sup F ( '^5"^'^;^' < e-^"^ < e-^" 

xj^o V \\9fx\\ ) 

Proof. As already mentioned the proof is an adaption of |DSS021 IBL85| proposition VI. 2. 2. for the 
case of PR,^ to our setting. In order to bound the probability 

we are looking for a function that upper bounds the characteristic function x- A possible choice is 
given by 



1 0< z< e-^"- 

2 - ze^" e^=" <z<2e- 
2e-^" < z<l 



(53) r„(x) = /„(|(^,y)|), /„(z): 

Computing the desired probability we get from equation ()52p 
(54) 

F [ ^fgZl^ < < E(r„(5,"x)) < I I FrMl^zidg) - I Tn{x)iyz{dx)\ + | J Tn{x)iyz{dx)\ 

using the definition of the expectation value and the triangle inequality. The first term can be 
bounded by a theorem about the convergence of the invariant measure (see V.2.5 in |BL85| ). which 
states that with the given assumptions there is a oq > such that for all < a < there are finite 
constants < Cq and < pa < 1 such that 

(55) IT^iTnigfx)) - ^,(r„)| < WT^ivMx)) - i^z{Tn)\L < CaPl ||r,|U 
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holds. Noting that x E PC^ can be written as a; = (r, ^fl^^r^ e^'^) with < r < 1 , using the 
definition of F, the mean value theorem for / and the unitary invariance of the scalar product and 
5 we find 

(56) |r(xi) - r(x2)| < II (xi ,y)| - I (X2 ,y)||e^" < In - rale^" < <5(xi,X2)e^" . 
In the last step we used that with the normal form of x G PC^ from above we have 

(57) 5{xi,X2f > {n,/l^ - r2^ll^,f > {n-r2f 
Equation (I56p implies now the following bound on the a-norm of for < a < 1 

(58) ||r,||„ = ||r,||^ + ?n„(r,) < ^/2e^" + 1 . 

A bound on the second term in equation ()54p can be achieved by regularity properties of the invariant 
measure (see corollary [Al8ll . Setting B = {x € PC^; |(x , y )| < 2e-^"} we get 



(59) I T,{x)u,{dx) < u,{{xeP(D';\{x,y)\ 

< 2e-""}) 

f \{x,y)f 
Jb I {x,y)f 

r II 11''^ 



{x,yW 

Using corollary I A. 181 the integral on the right hand side of this equation is finite for /3 small enough. 
Combining ()54p with the two bound ([58|) and ()59p we find that for some > 0, ng and < a < 
there are Cq ,K>0 and < /?q < 1 such that 

(60) F ( ^^"^'r^^ < e"""^ < 2^Ke-'(^'' + CaPa{V2e''' + 1) 

V \\9fx\\ ) 

For < e < log p = eq we can find no and (5 > such that the right hand side is bounded by e~^^ 
for all n < TiQ. □ 

By using proposition IA.14l iii and corollary lA. 201 there also exists a Eq > such that 

I (5?x,y)| > ||5^x|| e"^" > e^^^^^'^^)" 

for all < e < Eq with probability larger then 1 — e^" — 6°^" for all n larger then some no, we have 
also proven corollary IA.19I 

Our second aim was to prove that the Lyapunov exponent 7(2;) is Holder continuous on T for the 
cases considered in this paper. We transfer the z-dependency from our measure [i^ to the matrices 
(7, more precisely, we assume there is a measure [i and z-dependent matrices such that we may 
write 

' j^g^iMa) = [ fi9z)Kdgz) =-^{f{gz)) 



for all bounded measurable /. We still impose the following integrability assumption on the measure 

(61) 3C>0VzGT: E (^115^11'^) < 00, 

which also implies that IE (log | Iff^ 1 1 ) is finite. Matrix products of length n, where each factor is given 
by a matrix gz are denoted by gf. Another assumption concerns the dependency of the matrices g^ 
on the spectral parameter z €T, we assume that 

(62) ]S{\\gz>g~'\\)<C\z-z'\, z,z'gT, 

which is certainly satisfied by the matrices Tz{U) in remark [A. 121 
The Lyapunov exponent admits the following representation 

7(z) = M<^z), 
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with the function 

$,(x) = Eflog%^ 

and the unique invariant continuous measure v^, corresponding to /i^ , see lemma IA.14I The following 
proposition concerns some continuity results about $2 and 7(z) and is a summary of propositions 
V.4.8 and V.4.9 in [ULQO] . 

Proposition A. 21. Given that [ijLz) is strongly irreducible and non-compact together with assumption 
([62|) implies that 

i) the function $2 is continuous on PC^ x T and Holder continuous on T 

ii) the Lyapunov- exponent 7(2;) is continuous in z and the convergence 

j{z) = lim -E log-^^^^ 

n->oo n \ W^W 

is uniform in z and x. 

Proof. By lemma P\..14l <l>, is continuous on P<C^ for fixed z, hence, it remains to prove the Holder 
continuity of $2(2;) in z with constants independent of x. So, i) follows from 



|$2(S)-'J>.'(S)| <Eflog^^) <E(log||<7.'5.^^||) <C7| 

\ I l^z^l I / 



z-z'\ 



where the last inequality is just assumption (j62p . 

In order to prove ii) note that under our assumptions on fi^ it is guaranteed that the invariant 
measure corresponding to /iz is unique, cf. |CL901 IBL85| . This already implies that the Vz is 
weakly-continuous in z. If a sequence z„ G T is converging to A, then 

lim j{zn) = lim i/2„($(z„)) 

71— ^OO 71— S-OO 

= limiyz^{^z) + i^zA^z^-^z) 

n— ^-oo 

= M^z) 

= liz) 

Hence, 7(2) is continuous on T. 

If we define probability measures i^n,z,x = ^ Ylk=o l^^^x choose sequences z„ S T converging 
to z and Xn G PC^ converging to x we see that the probability measures Vn,zn,x„ converge weakly 
to the invariant measure v^. Clearly, we have 

hn{z,x) := -'E{\oggzX) = l'n,z,x{^z) , 

n 

and moreover hn{z,x) is continuous on T x PC^ by arguments similar to those used to prove i). 
Since T x PC^ is compact it is sufficient to prove that hn{zn,Xn) converges to 7(2). This follows 
from 



lim hn{Zn,Xn) = lim l^n,zn.x„{^zj 
71— >oo n— >oo 

= lim {l^n,z^,xA^z) + yn,z^,xS^Zr. " ^ z)) 
n— ^oo 

= Vz{^z) 
= 7(2) 

□ 

The Holder continuity of the Lyapunov exponent follows from the simple decomposition 
(63) 7(z) - 7(z') = - iyz')i<^z') + M'^z - <^z') 

if we prove that the statement is true for both parts of the sum. This is easy for the second part, 
since it was shown in proposition IA.21I that the map $2 is Holder continuous in z independent of 
the argument in PC^. So, we are left to prove the Holder continuity of the first term in 



36 ANDRE AHLBRECHT, VOLKHER B. SCHOLZ, ALBERT H. WERNER 

For that, let us first define the following distance on PC^, 

_^ \det{x,y)\ 

IfII lly|l 

where the determinant of two vectors x,y G is understood as the determinant of the two- 
dimensional matrix with rows x and y. The map defined in (j43p maps Ca to itself, and as a 
mapping defined on a Banach space it inherits some nice properties. We consider R^ as an element 
of B{Ca), the Banach algebra of bounded operators on with respect to the operator norm, which 
we denote by ||M||^ for M G B{Ca)- The following lemma provides all necessary tools for the proof 
of the Holder continuity of j{z). 

Lemma A. 22. For z € T define the operator Qz = Rz — ^z, then 
i) For G C with \9\> max(l,lim sup HQ" 1 1 ) resolvent of Rz exists and is given by 

-1 Nz , ^ 



Gz,9 := (Rz - Oiy = + 



1 ^ 0n+l ■ 
n=0 

ii) There exist positive constants Ca and pa < 1 such that \ \Qz\\^ < CaPa- Moreover, for 9 G C 
with \9\ > 1 we have the uniform resolvent bound 

iir II < 1 + 

iii) For e > denote by the circle of radius 1 + e around z = 0. Then we have the Cauchy 
identity 

1 



iV. + l = ^y Gz,ed9. 

We postpone the proof of this lemma to the end of this section. First, we prove the main result 
of this section, the Holder continuity of 7(2:). 

Lemma A. 23. The Lyapunov- exponent 7(2) is Holder- continuous on T, i.e. there exist positive 
constants a > and Fq, such that 

\l{z)-^{z')\<T^\z-z'\' . 

Proof. According to the discussion above it is sufficient to prove that the operator Nz is Holder 
continuous in z. We already know ^z £ C,a for arbitrary a > 0, hence we can apply lemma lA. 221 to 
see that 

A, - A,, = ^ y {Gz,e - Gz',e)d9 , 

Be 

hence, we may estimate 

(64) \Nz{f)-Nz'if)\ = \\Nz{f)-Nz'if)\\^<^ J \\Gz',e{Rz'-Rz)GzAf)\\j(^- 

Be 

The crucial point is to establish a Holder-estimate of the norm 

(65) \\{Rz' - Rz){h)\\^ = \\{Rz' - RzmW^ + m^{{Rz> - Rz){h)) , 

which we do for both terms separately. For /3 > and h ^ Cp the first term obeys the estimate 

\\{Rz' - Rz){h)\\^ < mp{h)snv^(5{gz'X,gzxf^ 

< m^(/i)E (^llc/^/ll^llg^ll^^ ■\<lei{gzix,gzx)f 

(66) < mp{h)^(\\gz>\f\\gz\f)-\z' -z\f' , 
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where we used ||(7z2;|| > Hg^H"^ for x E PC^ and Idet^^l = 1 at the first inequahty and (j62p 
in the last step. The integrability assumption ()6ip assures that there exists /? > such that the 
expectation value in (|66p is finite, i.e. there exists a constant < oo such that \\{Rz' — Rz){h)\\^ < 
C^||/i||^|z' — z\^ for positive but sufficiently small /3. With such a /? at hand we choose a = /3/2. 
Spelled out, the second term in (|65p reads 

(67) ma[{Rz' - Rz){h)) = sup — 

^¥=y 6{x,y)2 

(68) |(i?,. - Rz)ih){x) - {Rz' - Rz){hm\ < A„||/i||2„min(|z' - 5 (x, y)^-) 
and together with the estimate min(C^/?/, y) < C this implies 

ma((/?.' - Rz){h)) < AM2a\^' - • 
The assertion follows from lemma lA. 221 ii) and ^q, sup ||<I>2||2q, < 00. □ 

Proof of lemma \A . 2^ The operators Nz and Qz are both bounded and satisfy QzNz = NzQz = 
and N'^ = Nz- Hence, it follows from the triangle inequality for the operator norm and the 
completeness of B{Ca) that the series 

^ {Nz + QzT Nz , A 



+ 



n=0 n=0 

converges to a bounded operator for \9\ > max(l,limsup ||Q"|| ), which proves i). 

Property ii) follows from a general theorem concerning cocyles, see proposition IV. 3. 15 in |CL90) . 
if we prove the existence of an integer N such that 

^«n\ f ^ ^{9zX,gzy) Nfn n ^ ^ 

(69) sup / log _ /X {dgz)<Q. 

Again, denoting n-fold products of matrices Qz by gf, we get the following relation for the distance 
6 

, ^{gz^^dzV) 1 II n-ii 1 II n-ii 

log = -log 5-^3; - log 5f^y 

d{x,y) 

which follows from |det((72)| = 1 and ||x|| = \\y\\ = 1. Bv IA.211 the two terms on the right hand 
side converge uniformly with respect to z € T and x, y G PC^ to 7(^2), and moreover 7(2;) is also 
continuous bv IA.211 Hence, there exists an integer N such that almost surely 

1 K9zX,gzy) ^ „ 
sup log < , 

z&J,xj^y 0{x,y) 

which proves ()69p and hence ii) is verified. 

Now, the Cauchy identity iii) follows immediately from i) and ii). □ 

Appendix B. Measures with nonempty interior 

This section deals with properties of transfer matrices induced by probability measures on U{2). 
As explained in section a measure fi induces a family of measures fiz on Und- Our aim is to 
analyze the support of the measures jiz for measures /x with nonempty interior and prove that the 
group {^z) generated by //^ is non-compact and admits no reducible subgroup of finite index. The 
requirements of lemma [A. 141 namelv that there is a > with IE (|a|~'') < 00, are fulfilled in this 
case, see remarks IA.12I and IB. 31 The next proposition proves that (/z^) to be non-compact if the 
support of /i has nonempty interior. 

Proposition B.l. Let fi be a measure on U{2) such that supp(/x) has nonempty interior in the 
standard topology ofU{2), i.e. there exists U G supp(/Li) and an open set Ojj C supp(/.i) with U G Ofj. 
Then {^z) is non-compact for all z G C\{0}. 
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Proof. According to lemma there exist r G ]R,+ and a, /3,-y G [0, 2tt) such that the transfer matrix 
corresponding to U € Ojj and z € C\{0} can be written as 

Since and its inverse are continuous maps, we have that Tz{Ou) is an open subset of supp(^z). 
Hence, there exists r' ^ r such that T{r' ,a, /3,'y) G supp(/i). A straightforward calculation shows 
that 

/ f( /^ 9ir,r') \ 

(70) T(r, a, /3, 7) •T(r', a, /3,7)-i = «\ I"' 

with 



/(r, r') = +r2)(l +r'2) - rr' 



g{r, r') = r\/ 1 + r'^ — r'\/ 1 + . 
Hence, the eigenvalues of this matrix are given by 

A± = f{r,r') ±g{r,r') . 

Since the determinant of (f70|) is of modulus one these eigenvalues have to obey |A+A_| = 1. Therefore, 
either |A-|-| = 1 or there is an eigenvalue of (|70p with modulus strictly larger than 1, in which case 
{fiz) is clearly non-compact. 

Hence, it remains to prove that |A-|-| 7^ 1. This follows from the equation 

P+l - l-^-ll = 2 • m.m{\f{r,r')\,\g{r,r')\} , 

and the observation that f{r,r') = is not possible and g{r,r') = leads to r = r' . □ 

Next we show that non emptiness of 1^1 z implies strong irreducibility of (^^). 

Proposition B.2. Let fi be a measure on the group of unitaries U (2) such that the interior of supp{p) 
is nonempty. ThenMz G C\{0} we have that all subgroups of finite index in {fj,z)are irreducible. 

Proof. We choose an open set O G supp(//) and consider the corresponding set Tz{0) C GL(C,2). 
Since and its inverse are continuous Tz{0) is also open in the standard topology of GL(C, 2). For 
M G Tz{0) the set Tz{0) ■ is an open neighborhood of 1, hence there exists e > such that 

{1 + 6U : 0<6<e,U e U{2)} C r^(0) • M"^ . 

Therefore, for arbitrary normalized vector v we have 

Be{v) C Tz{0) ■ ■ V , 

which implies that the set Tz{0) • w has nonempty interior for all normalized tf; G C^. This already 
excludes strong irreducibility, since the interior of a finite union of subspaces is empty. 

□ 

Remark B.3. For measures ^ having a component which is absolutely continuous with respect to 
the Haar measure H the support supp(//) has nonempty interior. Moreover, for all C, > Q we have 

(i«r^) = (i«r^) + ^i^s (kr^ 

where fic resp. denote the absolutely continuous resp. singular part of the measure fx. The first 
term can be estimated by the expectation value with respect to the Haar measure 

E,. (|ar^) < ||/..|LE^ (lar^) = \W\L [ , 

which is finite for all < ( < 1. The second term is finite if the singular measure consists of 
finitely many point^. 



'Note, that we can assume that each of these points fulfiUs \a\ > 0, see section WM 
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Appendix C. Discrete measures 

In this section we give a proof of corollary 13.41 that establishes the dynamical localization of a 
disordered quantum walk, where the coin operations are drawn from the set 



with a, 6 € C fixed, but arbitrary, such that < |ap + \b\'^ = 1 and \a\ < \3b\. The actual probability 
distribution according to which the coins are drawn is irrelevant for our purposes. Our aim is to 
prove that the group generated by the corresponding transfer matrices satisfies the requirements of 
theorem 13. II for all 6q G T, i.e. we have to prove that {hoq) is C-integrable, non-compact and possesses 
no reducible subgroup of finite index. The first requirement holds easily for a ^ 0. We proceed by 
proving that (fi^) is non-compact for all z € C\{0}. 

Proposition C.l. Let ^ he a measure inU{2) such that supp(/i) = {H, X}, then (fi^) is non-compact 
for allz£ C\{0}. 

Proof. The transfer matrices corresponding to H and X are given by 

T„ = rAH) = ( -f ] a„d T. = rAX) = i ( _\ 



Similarly to the proof of proposition IB. II we consider the product 

1 _ 1 / b-V2 (1 - V2b) 



1 



^^■^^ " a V {V2b-l)z V2- 
and determine the modulus of its eigenvalues, which is given by 

,^ , \i3b±J\a\^-\3b\^\ 



\a\ 



Under the assumption \a\ < \3b\ there exists an eigenvalue X± with modulus strictly larger than 
one. □ 

In order to prove that a disordered quantum walk according to corollary 13.41 exhibits dynamical 
localization we have to prove that (nz) is strongly irreducible for z G C\{0}. This can be shown by 
using proposition II. 4. 3. in |BL85| . The original proof in |BL85| is given for GL(1R,, 2), therefore and 
for the convenience of the reader we repeat this proof for GL(C, 2). 

Proposition C.2. Let fi be a measure on GL((D,2) such that {fj,) is non-compact and \ detM| = 1 
for all X G supp(;u). Then (/i) is strongly irreducible if for any v in the projective space PC^ the set 
{X V : X £ {fi)} contains more than two elements. 

Proof. Suppose there are elements vi . . .Vn G PC^ such that (/i) vi L) . . . L) Vn = vi L) . . . U Vn- 
Each g £ (fj,) induces a permutation iTg on the set {vi, . . . Vn} and the map (7 — )• Tr^ is a group 
homomorphism. Denoting the kernel of this homomorphism by K we can construct the quotient 
group G/K which is isomorphic to the symmetric group Sn, hence G/K \s finite. Now, every g £ K 
acts on the corresponding vectors Vi G as Mvi = XiVi. Assuming n > 2 we can choose three 
elements vi, V2 and U3 and find coefficients a, /3 G C with a 7^ 7^ /3 such that = avi + /3v2- 
Applying g £ K yields the equation 

Xsavi + A3/5f2 = Xiavi + Ai/3u2 , 

which implies Ai = A2 = A3 = A. Hence, g = XI for all g £ K and since | det g\ = 1 this implies 
A = e^"^ with G R. This is a contradiction to the non-compactness of (^u) which excludes that both 
K and G/K axe compact. □ 

It remains to prove that the transfer matrices corresponding to }{ and X generate more than two 
elements in PC^ when applied to an arbitrary element v £ PC^. We prove a slightly weaker, but 
sufficient result. In order to apply lemma II. 141 we might show that (/i^) is non-compact and strongly 
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irreducible for all but finitely many z G T, hence it suffices to prove strong irreducibility for almost 
ah z £ C 

Proposition C.3. Let v £ PC^, z G C2\{0} and {fiz)the group generated by Th and Tx- Then 
{fJ,z)v contains more than two elements for all but at most two z G C\{0}. 

Proof. The elements v G P<C^ are given by the vectors 

In order to check whether there are three different elements in {fiz)vr,(j> we first observe that this 
is the case if v,.^^ is not one of the two eigenvectors of Th ■ T^^ because then Vr,,j>, Th ■ T^^Vr,^ 
and {Th ■ T^^ Vr^^ are linearly independent. The same argument applies for the eigenvectors of 
T^^ ■ Tx, hence we have to analyze the cases when some eigenvectors of Th ■ T^^ and T^^ ■ Tx 
coincide. The condition that one eigenvector of Th ■ T^^ coincides with one eigenvector of T^^ ■ Tx 
gives us a relation depending on b and z. There are four such relations corresponding to the four 
possible choices of signs in 

d\b-V2±ij _ 2^b+V2±i'y 
V2b - 1 ~ ^ V26 + 1 

with 7 = y^J6p — |ap G R. Note that b = ibl/-v/2 is excluded by the condition \a\ < 3b. There 
are at most two solutions for z to this polynomial equation, hence there are at most two z such that 
(fiz) is not strongly irreducible. □ 



Appendix D. Absence of continuous spectrum 

For the sake of completeness we include in this chapter the proof of a unitary and discrete time 
version of the RAGE theorem |Gol85| . For this purpose we follow along the lines of a proof for the 
Hamiltonian case, where the time evolution is induced by a self adjoint operator |Tes091 lKir07| . The 
goal of this section is therefore to connect the spectral properties of a realization of a family the 
disordered walk operators {W^j} with the dynamical behavior of vectors of the Hilbert space. 

In order to do so, we are interested in the properties of the measure pZ'^{t) of the time evolution 
of a realization VF^ of a disordered quantum walk quantum walk 

p-'yit) = {6y^cP, W% ®^)= f . 

We can at first prove the following version of Wiener's theorem: 

Theorem D.l. Let p be a complex Borel measure on the unit circle T and define 

fi{t) = [ 6*p{de) 
Jt 

then the time average of fi{t) has the limit 

t=o AeT 

Proof. Starting from the definition of fi{t), using linearity of the integral and the that ^ G T we get 

1 ^ 
lim y = 

t=o 

= ^lini^ r r (^J2^'^''~'A me'Widiey)). 
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Since the geometric series within the parenthesis is bounded by T and converges point wise to the in- 
dicator function X{o}(^) can interchange the hmit with the integration by dominated convergence 
and arrive at 



— 7r J — 7r 



r K{y})n*{dy)= J2 IMKDI'- 

J —n _r \ 



□ 

For a given unitary W we can decompose the Hilbert space into three orthogonal subspaces 
T^ac, T~Lsc and T-Lpp each containing the vectors (j) for which the spectral measure is absolutely 
continuous, singular continuous or pure point, respectively and these subspaces are left invariant by 
W . In addition define H.^ = T~Lac ® T~Lsc- A connection between this decomposition and the time 
evolution of the quantum walk is given by the following theorem: 

Theorem D.2. Let W he a unitary operator and G a compact operator then for ipc G T~ic 

T 



lim — ^VIIGVFVcII^ = 
t=o 



holds. 



Proof. Choose a vector '0c G T~ic- Since W leaves He invariant, by {(f), Wip) = {P}{^(j) ,Wip) we see 
that if is continuous so is p^,^. For the compact operator G there is a sequence of finite rank 
operators G„ = X^^^q "^^ i^l' ' • )^ converging to it. By the triangle inequality we only have to check 
the single rank one summands so we end up with 



T 



t=0 t=0 
which follows from Wiener's theorem. □ 

Equipped with this result, we can prove the discrete time version of the Rage theorem: 

Theorem D.3 (RAGE). Let W be a unitary operator and Gn a sequence of compact operators 
converging strongly to the identity. Then we have for He and Hpp 



1 

= {0 e ^ ; lim lim — — V \\GnW'il;\ |' = 0} 

n-5-oo T-s-oo i + 1 ^^-^ 
t=0 

npp = {^€'H; lim sup 1 1 (1 - G„)W^*0| |^ = 0} . 



Proof. We start with the continuous case, which by theorem ID. 21 holds for every ip G T-L^. Decom- 
posing now an arbitrary (f) ^T-L into tpc G T-Lc and tpp G T-Lpp we can infer that we have to find a lower 
bound on 



T 

lim lim y"\\GnW^^„\? + \\GnW^^c\? -2\{GnW^^„,GnW^^c 



i=0 

Since the norm of the G„ is uniformly bounded, because they converge to the identity strongly, 
and by theorem ID. 21 the last two summands tend to zero. So we have to show, that HGnVF^-^pH is 
bounded away from zero for n large enough. Instead we prove that 

lim sup 11(1- Gn)I^Vp| I = 



n— >oo 



t>0 
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Being an element of Tipp tpp can be decomposed into eigenvectors •i/'fc of the unitary operator W. 
Inserting this decomposition for ^l>p we can upper bound the norm by 

TV oo 

hm sup V|afce-^^'=| 11(1- G„)V'fc|| + 111- Gnll V |aA,.e"^^'^| HV-fcll 

- k=l k=N+l 

The first sum goes to zero by strong convergence of the Gn and the second goes to zero if we make 
large enough and using the fact that a strong convergent sequence of operators is bounded. At 

the same time this proves the second claim of the rage theorem for S Hpp- 

If we now again decompose an arbitrary vector ■!/; G in its components in He and Tipp as in the 

first case we are left to prove that supj>o ||(1 ~ G'„)W^*'0c|| stays strictly larger than zero for all n. 

Assuming the contrary we find 

1 ^ 

||H^Vc||< lim ——y\\{l-Gn)W'i^M + \\GnW% 



T^oo T + 1 ^ 
t=0 



<sup||(l-G„)I^Vc|| ^^0, 
t 

by strong convergence of the {Gn}- This contradiction then concludes the proof of the rage theorem. 

□ 

Finally we can prove the absence of continuous spectrum of a walk operator. 

Lemma D.4. Let Wi^ be a realization of a family of walk operators {W^;} that exhibits dynamical 
localization, e.g. the assumptions of theorem \3. 1\ can be verified for the complete unit circle T. Then 
with probability one W^j does not have continuous spectrum. 

Remark D.5. Lemma \D.4\ o-^so holds for the open arc of the unit circle L^ under the assumptions 
of theorem \3.1\ The proof is essentially the same as the one given for lemma D.4 bellow, with xi^s) 
inserted before all W^. 

Proof. In order to prove the statement we show, that for a given realization W^j of the family of walk 
operators {W^}, the subspace Tiduj) contains with probability one only the zero vector. Therefore 
we look at the projection of an arbitrary ip £ Ti onto T-Lc{uj). Choosing the projectors G„ = 
I]|«|<n Si -P-Si^ei, where P5,®e, denotes the projector onto the state 6i (g) and by HV'H = 
we get 



T 



(71) 



r + 1 ^ 

t=0 



-Gn)W^Pn^^^)^P\\' + \\GnW'^Pn^^^)i;\f 



T + 1^ 

t=o 



T 

^ "'1 - Gn)Wii;\\' - 11(1 - Gn)W^Pn^^^^^i;\f 



Note that the G„ are a sequence of operators of finite range that converge strongly to the identity. 
Taking the limit with respect to first T and then n of eq. (j7ip we see that the sums over the second 
and the third term on the right hand side tend to zero by theorem ID. 31 

To make the connection to theorem 13.11 we take the expectancy with respect to w on both sides 
and get with Fatou's lemma 




E (I li'H.H V'l r j = E ( jim ^lim J] 1 1 (1 - G„)I^^^ 

(72) 



<liniinfE hm ^ 1 1 (1 - G^M^^^I T 
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In order to proof that Tiduj) only contains the zero vector almost surely, we have to check that the 
right hand side vanishes for every ip € Ti. We do this by checking it for the total set {6x ® ei,x G 
= 1,2} of our one dimensional lattice. 
Inserting the definition of the G„ into equation (|72p we find for ijj = 6x <^ ei that 

liminf E ( hm V 1 1 (1 - G„)T^^,5, e^j 



t=o 

2 



< lim inf E > > sup 

\|i|=n+ij = i - 
oo 2 / 

< liminf > > E sup \ {5i ej ,W 5x ^ ei)\' 

\l\=n+lj=l 



\l\=n+lj= 
oo 2 

< liminf y y Ci(e-^2l'+"l'+e~^2l'-"l') 

l=n+l j=l 

holds because of the dynamical localization of the disordered quantum walk, i.e. lemma [4.251 For n 
large enough / is larger than |x|, so we can forget about the absolute values in the exponents and the 
sum takes the simple form of a harmonic series and we find for the projection of 5x <^ Cj onto Tiduj) 

oo 

E(||P^^(^)5,®e,||^) < liminf 2C7i(e^2^'+e-^^^') V e'^^'' 

\ ' ' ^ ' ' ' / n— >oo ^ — ' 

l=n+l 

The fact that e''^^^^ < 1 ensures the convergence of the harmonic series and by throwing away the first 
n terms and making n large, the expression tends to zero as required. Positivity of | \ P-}{^(uj)^x (i^ 
then already implies Py^^(^^^6x (8) = for all x and i almost surely and totality of the set of all 
5x Cj then completes the proof. □ 
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